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Abstract 

We classify the parity- and time-reversal-violating operators involving quark and gluon 
fields that have effective dimension six: the quark electric dipole moment, the quark and 
gluon chromo-electric dipole moments, and four four-quark operators. We construct the 
effective chiral Lagrangian with hadronic and electromagnetic interactions that originate 
from them and serves as the basis for calculations of low-energy observables. The form of 
the effective interactions depends on the chiral properties of these operators. We develop 
a power-counting scheme and calculate within this scheme, as an example, the parity- and 
time-reversal- violating pion-nucleon form factor. We also discuss the electric dipole moments 
of the nucleon and light nuclei. 



1 Introduction 



The Standard Model (SM) of particle physics contains, in its minimal version, two sources of 
time-reversal (T), or, equivalently, CP violation. In the electroweak sector, the phase in the 
quark mixing matrix pQ is associated with the Jarlskog parameter Jcp — 3 x 10~ 5 [2]. The 
strong sector contains the QCD vacuum angle 9 [3j, but the experimental upper limit on the 
neutron electric dipole moment (EDM) shows that 9 is unnaturally small, 9 < 10 -10 CP 
violation within the SM is believed to be insufficient for a successful baryogenesis scenario [5], 
and therefore new sources of CP violation are expected in order to explain the cosmological 
matter-antimatter asymmetry. This is not a surprise since the SM is likely but the dimension- 
four part of an effective field theory (EFT) that contains higher-dimensional operators, some of 
which will violate CP. 

Powerful probes of such CP violation beyond the SM are EDMs of nucleons, nuclei, atoms, 
and molecules [61 [7] , which violate both parity and time reversal (-PJT) . Since the SM predictions 
from the quark mixing matrix j8] are orders of magnitude away from current experimental 
limits, a finite EDM in upcoming experiments would be an unambiguous sign of new physics. In 
addition to impressive improvements [9] on the time-honored EDM experiments with neutrons 
and neutral atoms, in particular 199 Hg, which have resulted in very precise limits [10} lllj. novel 
ideas exist for the measurement of EDMs of charged particles, such as the proton, deuteron 
and helion, in storage rings [12]. An important question that comes up is whether, when future 
experiments measure nonzero EDMs, we will be able to pinpoint the microscopic source of P 
and T violation. 

EDMs of strongly interacting particles arise from the higher-dimensional operators at the 
quark-gluon level. These non-renormalizable operators might have their origin in a renormaliz- 
able theory at a higher-energy scale, such as, for example, supersymmetric (SUSY) extensions of 
the SM. At the SM scale, the most important higher-dimensional f^f operators should be those 
of dimension six, as we are not concerned here with CP violation in the dimension-five leptonic 
operator [13J that gives rise to neutrino masses and mixings. From symmetry considerations it 
is found p~H HS1 HU H7] that the following flavor-diagonal f*p operators appear at an effective 
dimension six: the quark electric dipole moment (qEDM) [18], which couples quarks and pho- 
tons; the quark chromo-electric dipole moment (qCEDM) [19j, which couples quarks and gluons; 
the Weinberg operator [20] , which couples three gluons and gives rise to a gluon chromo-electric 
dipole moment (gCEDM) [2T] . and four four-quark operators [TBI 1^2"]- 

Since it is not feasible to calculate hadronic and nuclear properties directly from a Lagrangian 
at the quark and gluon level, we use chiral EFT [23] — a generalization to more than one nucleon 
of chiral perturbation theory (%PT) [23] — to translate microscopic operators into operators that 
include nucleons, pions, and photons. (For reviews, see Refs. [25 1 126 1 [27].) After the translation, 
we are able to calculate hadronic properties directly from the effective Lagrangian. For the 
dimension- four 9 term, this method was first employed in Ref. [4], and later extended in the 
context of SU{2) x SU{2) [281 E3 M, ED E2J E31 Ell E3 E3 ED EE] and of SU(3) x SU(3) [391 ED] 

In this paper we extend the method further to include dimension-six operators in the frame- 
work of SU(2) x SU{2) xPT. (Generalization to SU(S) x SU(3) is straightforward.) Different 
sources at the quark-gluon level produce different hadronic interactions. The effective chiral 
Lagrangian includes not only interactions that stem from spontaneous chiral symmetry breaking 
and are therefore chiral invariant, but also interactions that break chiral symmetry in the same 
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way as chiral-symmetry-breaking operators at the QCD level. Since the dimension-six opera- 
tors break chiral symmetry differently from each other and from the 9 term, they will generate 
different low-energy hadronic interactions. Given enough observables it should be thus possible 
to separate the various f*p sources. 

In addition to constructing the Lagrangian, we need to organize in leading order (LO), next- 
to- leading order (NLO), etc. the various effective ff operators that appear. This is done 
according to the estimated size of their contributions to observables. In order to get a consistent, 
manifest power counting we work in the heavy-baryon framework |41j wherein the nucleon mass 
has been eliminated from the nucleon propagator. This framework has a transparent power 
counting and greatly simplifies loop calculations, but there are some complications when one 
goes to subleading orders in the Lagrangian. These problems can be solved by demanding that 
the Lagrangian obeys reparametrization invariancc (RPI) [42 i. This puts constraints on certain 
coefficients of operators, which we construct up to NNLO. 

In Refs. [H2 [331 EH EHl [36] we have used some of the effective interactions to calculate 
the EDMs of the nucleons and the lightest nuclei. It was found that in LO they depend on 
six low-energy constants (LECs), which for each fundamental p*p source have different relative 
sizes. The idea is that, from EDM measurements, the LECs can be inferred and the dominant 
fundamental source identified. For example, if the deuteron EDM is significantly larger than the 
sum of the nucleon EDMs, this points towards new physics in form of a qCEDM |34| [36] or, as 
we will demonstrate here, a particular isospin-breaking pf four-quark operator. If the deuteron 
EDM is well approximated by the nucleon EDMs, but the helion (triton) EDM is far away from 
the neutron (proton) EDM, it is a hint for the SM 9 term. The deuteron magnetic quadrupole 
moment (MQM), if experimentally accessible, could play an important role as well \64\ I37| . 

In this article, we construct the full Lagrangian necessary to perform these and other low- 
energy calculations. We follow the approach outlined in Refs. [30, [32], where the xPT 
Lagrangian originating from the 9 term was derived, and the higher-dimensional interactions 
were examined. We identify, for each source, the size of the six^^T interactions relevant for light 
nuclear EDMs. Apart from that, we construct also all other operators that appear at the same 
order. These operators could play a role in the calculation of other observables that violate P 
and T such as the EDMs of heavier nuclei or of ions and atoms, form factors and scattering 
observables could depend on these operators as well. An important ingredient in calculating 
hadronic and nuclear pf observables is the pf pion-nucleon form factor. We investigate this 
form factor here. Finally, we extend our earlier light-nuclear EDM calculations to a four-quark 
operator that has recently been shown to arise below the SM scale from weak-boson exchange 

E2]. 

Our paper is organized as follows. In Secs.[2]and[3]we discuss the QCD Lagrangian, the 9 term, 
and the possible higher-dimensional pf operators. In Sec. [4] we briefly discuss SU(2) x SU(2) 
^PT, and how to incorporate the pf operators in this framework. The bulk of the article is 
Sees. [5j [6j and [8j where we construct the xPT Lagrangian up to NNLO including the removal 
of pion tadpoles that proliferate when pf appears together with isospin violation. In Sec. [7] we 
use the constructed Lagrangian to calculate the pf pion-nucleon form factor (PNFF), and in 
Sec. [9] we give the electric dipole form factor (EDFF) of the nucleon in case of the specific pf 
four-quark operator of Ref. [22] . In Sec. 10 the extra elements that arise in nuclear systems are 
addressed. We discuss our results and conclude in Sec. dU 
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2 The underlying quark-gluon Lagrangian 



Before we discuss higher-dimensional operators, we set our notation by recalling the main SM 
ingredients. The SM Lagrangian is completely determined by gauge symmetry with gauge group 
SU C (3) x SUl(2) x by the matter content, and by the requirement of renormalizability. 

In its minimal form, which we assume in this paper, the matter content consists of three 
generations of leptons and quarks, and one scalar doublet. Quark and lepton fields carry a 
generation index r = 1, 2, 3, which we will often leave implicit, running over the three generations 
of up-type quarks u = (u, c, t), down-type quarks d = (d, s, b), charged leptons e = (e, fi, r) and 
neutrinos v = (u e , u^, v T ). The left-handed fermions are doublets of SUl(2), 



while the right-handed fields ur, dR and e# are singlets. The field ip denotes an SUl(2) doublet of 
scalar fields ip J , J = 1,2. For convenience we define (p 1 = e IJ ip J * , where e IJ is the antisymmetric 
tensor in two dimensions (e 12 = +1). Left and right-handed quarks are in the fundamental 
representation of SU C (3), while leptons and scalars are singlets. We sometimes group the right- 
handed up- and down-type quarks in a doublet qR, and define quark doublets q = qL + Qr- The 
hypercharge assignments under the group Uy(l) are 1/6, 2/3, —1/3, —1/2, —1, and 1/2 for q^, 
ur, d-R, II, e R , and ip, respectively. 

We denote the gauge bosons associated with the gauge groups SU C (3), SUl(2), and Uy{l) by, 
respectively, W^, and B^, with a = 1, ... ,8 and i = 1,2, 3. Gauge invariance is most easily 
imposed by employing covariant combinations of the gauge and matter fields. The covariant 
derivative of matter fields is 

D» = d,- i| Gl\ a - i 9 - Wy - ig'YB^ (2) 

where g s , g, and g' are the SU C (3), SUl(2), and Uy(l) coupling constants; and X a /2 and r*/2 
are SU(3) and SU(2) generators, in the representation of the field on which the derivative acts. 
For example, for left-handed quarks A a and t 1 are, respectively, the Gell-Mann color and Pauli 
isospin matrices. The field strengths are 

Gfu, = d^G® — d^G^ — g s f abc G b ^G c u , (3) 
W; v = d»Wi-d v Wl-g^ k Wrwt (4) 
= dpBv-dvBn, (5) 

with f abc and e ijfc denoting the SU(3) and SU{2) structure constants. 

The SM Lagrangian is expressed in terms of all possible dimension-four gauge-invariant op- 
erators: 

Csm = -\(G%G a ^ + W; u W i ^ + B^B^) 

+q L ip qi + u R ift u R + d R ip d R + T L ip l L + e R ip e R + D^D^ip 

-q L Y u pu R - q L Y d ^d R - l L Y e <pe R + H.c. 
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where e pva ^ is the totally antisymmetric symbol in four dimensions (e 0123 = 1). 

The first line of Eq. ^ contains the kinetic terms and self- interactions of the SU C (3), SUl(2), 
and gauge bosons. The second line contains the kinetic energy and the gauge couplings 

of fermions and scalars. These couplings are completely determined by gauge invariance. The 
terms of the third line form the scalar potential. With the parameter A > for stability and the 
parameter fj 2 > 0, the scalar field acquires a vacuum expectation value v = y / ^ 2 /A; the Higgs 
boson h(x) represents fluctuations around this vacuum, 



v 



v 



where U(x) is an SU(2) matrix, which encodes the three Goldstone bosons. The Goldstone 
bosons are not physical degrees of freedom, so that with a particular choice of gauge, the 
unitarity gauge, U (x) can be set to one. In this gauge, the Goldstone bosons are "eaten" by the 
longitudinal polarizations of the massive vector bosons. We will often use this gauge to discuss 
the structure of dimension-six operators. After electroweak symmetry breaking, the scalar field 
kinetic energy provides a mass term for the weak gauge bosons. It is convenient to express the 
fields W$ and in terms of the physical photon and Z-boson fields, 

W* = caaOwZp + innOwA-p, (8) 
Bp = cos 8 W - sm9 w Z fu (9) 

where the weak mixing angle 9\y is given, together with the proton charge e > 0, by the couplings 
g and g' via 

9=--^- (10) 



sin^' cos 9\y 

The charged I^-boson fields are defined as 

W± = -^(W M 1 T ^ 2 ). (11) 

The next dimension- four operators one can write are the Yukawa couplings of the fermions to 
the scalar boson in the fourth line of Eq. ([6]), via matrices Y u,d,e . After electroweak symmetry 
breaking, they generate the quark and lepton masses. By means of unitary transformations on 
the quark and lepton fields, it is always possible to make the fermion mass matrices diagonal and 
real, up to a common phase. For leptons, these transformations do not leave any trace, while for 
quarks the price to pay for the diagonalization of the mass matrix is that the interaction of the 
charged W bosons with the quarks is no longer flavor diagonal, an effect that can be obtained 
by replacing d r L in Eq. ^ by V rs d s L , where V rs is the unitary Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. For three generations of quarks, the CKM matrix has one complex phase pQ, 
which is responsible for the observed CP violation in the kaon and 5-meson systems. However, 
the contribution of the CKM phase to nuclear EDMs is orders of magnitude smaller than the 
current experimental sensitivity, and we will neglect it in the rest of the paper. The second 
J°JT parameter in the SM Lagrangian is the global phase of the quark mass matrices. It can 
be eliminated [H] by an axial rotation of all the quark fields q^ — > e ip qL, qn — > e~ lp qR. Such 
transformation is anomalous [15] , and its net effect is to shift the coefficients of the P- and T-odd 
gluon operator in the fifth line of Eq. ^ from 9 to 9 = 9 + njp, where nf is the number of 
quark flavors. Operators in that line are total derivatives, but, for non-Abelian gauge fields, they 
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contribute to the action through extended field configurations, instantons [5J. The contribution 
of instantons is proportional to exp(— 8n 2 /gf). For QCD instantons, the coupling constant g s is 
large at low energies, and the instanton contribution to the action is not negligible. On the other 
hand, the electroweak coupling g is small, and electroweak instantons are extremely suppressed, 
negligible for all practical purposes. We will neglect the electroweak theta term. 

The SM terms related to CP violation which we will be concerned with below can be sum- 
marized by 

A = -S^^G^Gap + (Mo + M 3 r 3 ) q R + H.c., (12) 

where Mo,3 are real, diagonal quark mass matrices. 

There are good experimental and theoretical reasons to believe that the SM is an EFT, at 
the electroweak scale Myy, of a more fundamental theory, and it is renormalizable in the more 
general sense of including all higher-dimensional operators allowed by the symmetries. In the 
form of Eq. (|6j) the SM does not include neutrino masses and mixings. These can be accounted 
for by introducing the only gauge-invariant dimension-five operator that can be written in terms 
of SM fields [13 . The lightness of neutrinos compared to other fermions is explained if the 
new-physics scale that suppresses this operator is of the order of 10 15 GeV. Any new physics 
between the electroweak and the GUT or Planck scales will manifest itself in higher-dimensional 
operators. A promising strategy to probe such operators is to look for processes where the SM 
contribution is extremely small. Examples are rare processes like the lepton flavor-changing 
process fj, — > ej or neutrinoless double beta decay. 

Our work focuses on EDMs, which signal CP violation in the flavor- diagonal sector and are 
insensitive to the phase of the CKM matrix. Since is very small, it is possible that higher- 
dimensional fp operators are competitive with the 9 term. These higher-dimensional operators 
might eventually be linked to an underlying, ultraviolet complete theory. We denote the 
scale characteristic of this theory by Mp. Well below the scale M-p we expect effects to 
be captured by the lowest-dimensional interactions among SM fields that respect the theory's 
gauge symmetry. In general, operators of dimension (4 + n) at the SM scale are suppressed by 
powers of M^ n . The next-to- lowest-dimension operators involving quark and gluon fields 
that can be added to the SM Lagrangian have effective dimension six. Their complete set was 
constructed in Refs. [141 1151 120 ] 116| and recently reviewed in Ref. |17j . 

We focus here on flavor-diagonal operators involving quarks and gauge bosons only. Fol- 
lowing Ref. [IT], we organize the dimension-six operators we need at the electroweak scale 
according to their field content: three vector bosons (A^ 3 ), two gauge bosons and two scalars 
(X 2 ip 2 ), two quarks, a scalar, and a vector boson (q 2 ipX), two quarks, two scalars, and a deriva- 
tive (q 2 ip 2 D), two quarks and three scalars (q 2 tp 3 ), and four quarks (q 4 )- The coefficients of 
these operators are all proportional to M^ 2 . 

At dimension six, the interactions of a quark with gauge bosons gain in complexity. With the 
aid of the scalar boson we can write 

^L^ v (f d \ a G% + T d B B, v + Yiyr'W^) ^d R + H.c. , (13) 



where, in the most general case, the couplings f u ' d and T B ' w , are 3x3 complex-valued matrices. 



After electroweak symmetry breaking, the operators in Eq. (13) can be expressed in terms of 
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dipole moment operators, 



C 



pv 



r + r 3 r 3 ) X a q R G% - -q L a^ (T + T 3 r 3 ) q R 



q L a^ u (Tzo + r Z3 r 3 ) q R Z, 



pv 



-l=d L a^T^u R W f 
where, at a renormalization scale /u ~ M^/, 



pv 



V2 



pv 







To,3 = 


i( 




1 " 


To,3 = 


2 . 




1 " 


Z0,3 = 


2 . 



± J cos W + (r^ =f r^J sin # w 
(rfr =f r&r) cos^ - (r£ ± r£) sm% 



(14) 

(15) 
(16) 
(17) 



The first two operators in Eq. (14) are the most interesting for low-energy applications. The 
imaginary parts of the diagonal entries of To i3 and ro >3 generate the quark electric and chromo- 
electric dipole moments (qEDM and qCEDM). Since these operators flip the chirality of the 
quark field, we assume these matrices to be proportional to the Yukawa couplings in the SM 
Lagrangian, and thus to the masses: 



0,3 



O 



r 



0,3 




(18) 



where 5$^ an d <5o,3 are dimensionless constants that parameterize any deviation from this as- 
sumption, and contain all the information on physics beyond the SM. The non-diagonal entries 
are also of considerable interest, since they produce flavor-changing neutral currents. In the 
CP-odd sector, for example, the uc entries were found [46: to be the less constrained dimension- 
six operators that contribute to the recently observed CP violation in charm decays [37] • Since 
the flavor-changing operators contribute to nuclear EDMs only via additional loops involving 
weak-boson exchange, for our current purposes we can neglect them. The next three operators in 
Eq. ( |14[ ) are weak dipole moments. Their contribution to the qEDM and qCEDM is suppressed 
by (g/47r) 2 . In the ". . ." in Eq. (14) we find couplings of the Higgs boson, which we can also 
neglect. 

The self-interaction of gauge bosons also exhibits new structures at dimension six. CP viola- 
tion can be found in 



C X 3 



dw rgbc pvafipa pb pep , jjk puaffriri xrfj wfc p 

g J c KJ al3^pp KJ u ~ g c c "o"" " -■ ' " 



pp' 



(19) 



The first operator is the Weinberg three-gluon operator [20] which can be interpreted as the gluon 
chromo-electric dipole moment (gCEDM) |21| . Similarly, the second operator is the T^-boson 
weak electric dipole moment. After electroweak symmetry breaking, this operator generates 
interactions containing at least two heavy gauge bosons [15], which can contribute to q(C)EDMs 
through loop corrections. Again such contributions are suppressed by (g/Ai:) 2 [15J. For our 
purposes we can neglect the weak dipole moment and focus on the gCEDM, 



w 



dw = 0\ 4vr- M2 



(20) 
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with w a dimensionless constant. 

Just as there is CP violation in multi-gluon operators, so there is in multi-quark operators. 
In agreement with Ref. [16], we find just two four-quark interactions, 

4,4 = Ex (q J L u R ) e JK (q£d R ) + S 8 {q J L X a u R ) e JK (qf \ a d R ) + H.c, (21) 

where the couplings Si 8 are four-index tensors in flavor space. They scale as 



O (V)^) , (22) 



c q 4 = 7 ( ImS i)im (.qqqijsq - qtq ■ qri^q) 



where a\ g are dimensionless constants. The operators in Eq. ( 21 ) are not affected by electroweak 
symmetry breaking. For later convenience, we rewrite Eq. (21) in terms of quark doublets 
q = qi + q R , and focus on the flavor-diagonal ff terms only, 

1 

4 

+ - A (Im£ 8 ) ml {q\ a qqi l5 \ a q - qr\ a q ■ qri l5 X a q) + ... (23) 

These four-quark operators, originating directly at the electroweak scale and constrained by 
exact SUl(2) gauge invariance, are chiral invariant. The remaining parts of Eq. (21) and other 
dimension-six four-quark operators listed in Ref. [T7] necessarily involve flavor-changing 
effects, and contribute to flavor-diagonal CP violation only at the loop level. We will neglect 
them in what follows. 

The remaining fp sources involve two or more scalar fields. The most important of them is 

C q2tp 2 D = u R Ei^d R $iDn<p + H.c, (24) 

where in general Hi is a complex 3x3 matrix in flavor space. Here we focus, again, on the 
flavor-diagonal parts only, given after electroweak symmetry breaking by 

2 

CqWD = 2^ W""« H iVdfl + W^d R Erfu R ) +... (25) 



Contrary to the operators in Eq. (13), the operators in Eq. (25) do not change chirality and we 
do not expect them to be proportional to the quark mass, so we parameterize 

Hi = f (47r) 2 -|a | • (26) 




At low energy, after we integrate out the W boson, the imaginary part of Hi contributes to 
four-quark operators, which are particularly interesting because they are not suppressed by 
the light quark masses [22]. As detailed in Sec. [3| ImHi generates four-quark operators of the 
same importance as the chiral-invariant four-quark operators that are generated directly at the 
electroweak scale, which we introduced in Eq. (21). 

Finally, there are other terms that are closely related to interactions in Eq. dm), 

£xV,<?V = ~ 32^2" \9b& G % G aP +9 Q'w W liV W aj3 + d' d 'b B^B^) 



327r^ ^ v v 



t » t 

W^Baf,?-^- - 2^ [q L Y' u ^u R + q L Y' d vd R ) , (27) 
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where the angles 9', 9' w b wb and the Yukawa couplings Y' u,d , which are symmetric, complex 
matrices in flavor space, scale as 




rj __ /n I I -\s>u,d 

w, b, wb 




Y ' U4 = tjo • (28) 



Equation (27) is reminiscent of Eq. Q. Indeed, if one rewrites 2^ip/v 2 = 1 + (2^ip/v 2 - 1), 
the pieces of the first three and last two operators that do not contain the Higgs boson can 
be absorbed in a redefinition of the couplings 9, 9 W: b, and Y u,d in Eq. Q. In some sense this 
aggravates the strong CP problem since even if the bare QCD vacuum angle is tuned to zero, 
we need to explain why the contributions induced by these higher-dimensional terms are small 



as well. Of the first four terms in Eq. (27), only the operator 6' b obtains after electroweak 
breaking a nontopological piece without a Higgs boson |15| . This piece connects three or more 
electroweak gauge bosons and contributes through loop diagrams to the quark electric and weak 
dipole moments, but such contributions are suppressed by (g/Air) 2 . The remainder of Eq. (27) 
then includes CP-odd operators with at least one Higgs boson: in the unitarity gauge, 

£*V,9V = - { 32^2 (9s°' G % Qa afi + 9 2 0' w W; U W^ + g' 2 9' b B^B a0 + gg'9' wb W^B a p) 

+V2v (l + -) (u L Y' u u R + d L Y' d d R ) \ - ( 1 + AV (29 ) 



v 1 V ' ' ' I v V 2v 



In the first line of Eq. (29) are CP-odd interactions of the Higgs to two gluons, two photons, 
or two weak bosons. The imaginary parts of the Yukawa couplings generate flavor- diagonal 
and flavor-changing Higgs-quark interactions. At low energy, where the Higgs boson is 
integrated out, these interactions manifest themselves in loop corrections to the qEDM, qCEDM 
and gCEDM. At tree level, PJT Higgs-quark interactions generate PJT four-quark operators, 
which, as discussed below, are relatively small. 

Of course, we could continue the procedure to higher-dimensional PJT sources. Those would 
include, for example, further four-quark operators [H] at dimension eight. We limit ourselves 
here to a systematic study of the effects of dimension-six sources to low-energy observables. 



3 Tree-level matching onto the QCD scale 



The PJT Lagrangian of dimension up to six that is relevant for the calculation of hadronic and 
nuclear EDMs is summarized in Eqs. (1121) , (14), (19), (23), (25), and (29). For low-energy 



applications, it is important to evolve the Lagrangian from the electroweak scale down to 
the typical hadronic scale [i ~ Mqcd — 1 GeV. In the process, one has to integrate out the effects 
of heavy SM particles |15|l49j; at the same time, one has to evaluate the running of the coupling 
constants and account for the possible mixing of the dimension-six operators |21l [50l I48j . A 
detailed account of the matching and evolution of the complete dimension-six P^T Lagrangian is 
beyond the scope of this work. Here we limit ourselves to tree-level matching, and we turn off 
(most of) the running of the couplings. 

At the QCD scale, the theory involves only light quarks, gluons and photons. Focusing on 
the first generation, 



C, 



cl 



[G%G a ^+F^F^)+q 



' + 



(30) 
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has a global, chiral SUl(2) x SUr(2) ~ 50(4) symmetry under independent SU(2) rotations of 
left- and right-handed quarks. Remaining dimension-four (and most of the higher-dimensional) 
interactions break chiral symmetry explicitly. Because all the breaking is due to relatively small 
quantities, it is useful to classify interactions according to their transformation pattern under 
50(4). 

For the construction of f*p electromagnetic operators in the chiral Lagrangian, we will need 
to consider the chiral properties of the PT electromagnetic couplings of the quarks in 

C e = -eArf^Qq, (31) 

where 

<?4 + ? (32) 

is the quark charge matrix. Introducing a chiral-invariant 1^ and an antisymmetric tensor 
through 



we rewrite Eq. (|31|) as 

-aJ 1 - 

2 M V 3 

At low energy, we can neglect the effects of the Higgs and of the heavy quarks, and Eq. (jl2j) 
is rewritten as 

U = -q L Mq R - q R Mq L - G%G a a p , (35) 

where 

with real parameters m u d, or alternatively 

m u + m d m d - m u 

m= , e = . (37) 

2 m u + m d 

For a xPT treatment, it is more convenient to eliminate the 9 term with an axial U(l) rotation 
on the quark fields, and move all CP violation to the quark mass term. After vacuum alignment 
|44j . the QCD 6 term becomes, in the notation of Ref. [32], 

£ 4 = mr(6)S4-efhr- 1 (9)P 3 -m*sm6r- 1 (8)P 4 , (38) 

where we introduced two 50(4) vectors 

*- F-($ t ), (39) 

a function r(9) that goes to 1 in the limit of small 9, 

f I + e 2 tan 2 ±9~\ 1/2 - 9x , s 

'«= W K 1 + 0(9)1 (40) 
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and the parameter 

m u m d l-e 2 _ 

= = m. (41) 

m u + m d 2 

From these ingredients one can derive the form of low-energy interactions originating in the 
QCD vacuum angle, as presented, for example, in Ref. [32] . 

We will now consider the dimension-six Lagrangian, which breaks chiral symmetry in new 



ways. The light-flavor components of the dipole operators in Eq. ( 14 ) match onto the light-quark 
EDM and CEDM, 

£ q( C)EDM = ~\ <?^"V (do + d^qF^ - \ qia^j 5 (d + d 3 r 3 ) X a qG%, (42) 
where, at tree level, 

d ,3 = Im(r ,3)n, do,3 = Im(f 0,3)11- (43) 



The non-diagonal components of the qEDM and qCEDM operators in Eq. ( 14 ) correct Eq. (|43j) 
at the loop level. The weak EDMs also generate loop corrections to the qEDM and qCEDM. At 
tree level, they generate four-quark operators, which, however, involve one power of the quark 
momentum, and are thus dimension-seven. The couplings of such operators are suppressed by 
an additional factor rriu^/M^ compared to the couplings of the dimension-six operators that 
we keep. 



The gCEDM in Eq. (19) and the four-quark operators in Eq. (23) match onto themselves. As 



pointed out in Ref. [22], at tree level, the only non-trivial result comes from the operator in Eq. 

W 



(|25|). The largest effect comes from a W exchange between light quarks, the suppression by M% 



from the W propagator being compensated by the factor g 2 v 2 from the vertices. The resulting 
-PJT four-quark operator, when evolved from M\y to Mqcd> induces another four-quark operator 
of similar form but with additional color structure Q The latter does not appear directly at 
the electroweak scale due to its gauge-symmetry-breaking properties. We find the following f^f 
operators, which because of their left-right mixing we abbreviate as FQLR: 

£ L R = i Im(~i)ii V ud (u R ^d R d L ^^u L - d R ^u R u L ^^d L ) 

+i Im(E 8 )n V ud {u R 'fX a d R d L l^ a u L - d R ^X a u R u L ^X a d L ) , (44) 

where V ud ~ 1 is a CKM element. Here, (Hg)n is not an independent coupling, but it depends 
on (Sx)n and on QCD renormalization-group factors, which are calculable. Thus (Hi)n and 
(Sg)n both depend on the same dimensionless parameter £. In terms of light-quark doublets, 



we can rewrite Eq. ( 44 ) as 



£ LR = i im(~i)n e 3ij qf'fq qr j j^ 5 q + i Im(H 8 ) n e 3ij qr l ^\ a q qT j -f^ 5 X a q. (45) 
By a Fierz rearrangement we can further rewrite the first term of this equation as 

- Im(Si)n [2 (qq qij 5 T 3 q - qr 3 q qij 5 q) - 3 (q\ a q qi^T 3 X a q - qr 3 X a q g«7sA a g)] . (46) 
o 

This operator was studied in Ref. [51] in the framework of left-right models. Although a large 
list of four-quark operators is presented there, only one combination (On — O12 + 6O21 — 6O22, in 



lr This was pointed out to us by W. Dekens. We thank him for useful discussions on this subject. 
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their notation) is generated at the electroweak scale. This combination is identical to Eq. (44). 
It is interesting to point out that we could have coupled the operator in Eq. (25) to the left- 



handed lepton current. The operator created this way causes /T in f3 decay through contributions 
to the triple correlation ~ D J- (p e xp v ). In Ref. |22j it is argued that, with current experimental 
accuracy, the best limit on the FQLR operator comes from EDM experiments. 



All other effects have extra suppression at low energies. The operators in Eq. ( 29 ) contain 
at least one Higgs field, and generate effects at either tree or loop level. For example, at tree 
level, the operators in the second line contribute to four-quark operators. But, since the Higgs 
couples to the light quark mass, the resulting f*p operators have a suppression of at least the 
small ratio of the light-quark mass over the electroweak scale. 

In summary, the (effectively) dimension-six Lagrangian is 



1 



q (d + g?3T 3 ) ia^jsq F f 



1 



2 

6 



q{d + d 3 n)ia> 11/ l5 \ a q G a u 



faocava.fi/-~ia /~tb s~icp 



+ 



ImHi 3ij 



4 

ImSi 



e *3q T ^ qqT 3 lffil5q + 



IUK 



q T i ^\ a qqT j ~i ll \ a >y b q 



(qqqij5q ~ qrq ■ qri^q) + (q\ a qqi^X 



q - qr\ a q ■ qrij 5 X a q) , (47) 



where all coupling constants have been redefined in order to absorb effects from operator mixing 
and renormalization-group running, and to drop the generation indices in the case of four-quark 
operators. As discussed above, one can imagine additional ff four-quark operators, but such 
operators are suppressed by weak gauge couplings in the typical combination (g/Air) 2 , small 
off-diagonal CKM elements, and/or powers of m U) d/My/. We therefore expect the operators in 
Eqs. (38) and (47) to give rise to the dominant P and T violation in hadronic and nuclear 



systems at low energy. 



The redefined constants in Eq. (47) scale as 



^0,3 



=1,: 




0\ 4vr 



O 



m2 t. 



(4tt)Vi,; 



iw 




0\ 4tt 



(48) 



in terms of dimensionless numbers 5o,3, <5o,3> w, £, and a\^. The sizes of these dimensionless 
parameters depend on the exact mechanisms of electroweak and P and T breaking, and on 
the running to the low energies where non-perturbative QCD effects take over. The minimal 
assumption is that #0,3, ^0,3, w, £, and ai t s are 0(1), 0(g s /Air) : (D((g s /4:ir) 3 ), 0(1), and 0(1), 
respectively. However, they can be much smaller or larger, depending on the parameters encod- 
ing _f^T beyond the Standard Model. In the SM itself, where = M]y, ^0,3, 5o,3 and w are 
suppressed not only by the Jarlskog parameter [2] Jcp ~ 3 x 10~ 5 , but also by additional powers 
of small gauge coupling constants and ratios of quark-to-VF masses |521 [7], In certain supersym- 
metric models with various simplifying universality assumptions of a soft-breaking sector with a 
common scale Mstjsy ; one has Mp = Mstjsy and the dimensionless parameters are of the size of 
the minimal assumption times a factor which is 0G3], roughly, Aqp = (g s /^) 2 sin (neglecting 
electroweak parameters), with eft a phase encoding T violation. If in the soft-breaking sfermion 
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mass matrices non-diagonal terms are allowed, enhancements of the type m^jra^ ~ 10 3 or even 
rrit/m u ~ 10 become possible, although they are usually associated with other, smaller phases 
[7]. In other models, the relative sizes of the dimensionless parameters could be different still. 

Below the hadronic scale Mqcdj the dimension-six sources generate further effective 
interactions, which break chiral symmetry in their own ways. Introducing the SO (4) singlets 

j eabc ^iivaB r~ia rib /—rep 

T m = 4 (99 9*759 - Qi-q ■ qrij 5 q) = -S ■ P, 

1$ = ±(q\\qi l5 \ a q-qT\ a q.qTi l5 \ a q), (49) 

the 50(4) vectors 

-iqa^rq \ „ T . _ 1 ( qa^rq 



"-a v r *~ ^=H4>™ (50) 



and 



rrr _ 1 / -iqa^ lb T\ a q \ v = l( q^ v r\ a q . 



and the symmetric 50(4) tensors 

y(i) _ 1 Z' Q t1 ^QQ tJ 7^Q ~ Qr^^qqr^^q -eP^qr^qqr^^q \ , , 

4\ -e lKL qT K YqqT L -f^q qTj^q ■ qrj^q - qrj^q ■ qrj^q J 



and 



y(8) _ 1 / qr t 'j' J -X a qqT^ fl X a q - qT t -/' J -/ 5 \ a qqT^^ 5 X a q -e jkl qT k ^\ a qqT l -f^ b \ a q 

A\ —£ lkl qT k ^\ a qqT l "f^<j\ a q qr^\ a q ■ qr^^q — qr^^\ a q ■ qrj^X' 1 



we summarize the Lagrangian as 



(53) 



£ 6 = -doVHdaWb-^+Wa+d^^^ (54) 

We already emphasized that the dimensionless coefficients <5o,3, <5o,3, w, £, and ui^ that we 
have introduced, and their relative sizes, strongly depend on the particular high-energy model, 
making it difficult to compare the relative contributions of different Tf^ sources to the same 
observable in a way that is independent of the details of the physics at the high-energy scale 
Mp. We therefore adopt the approach that we construct the low-energy Lagrangian for 
each dimension-six source separately. These separate Lagrangians can be used to calculate 
for each source their contribution to hadronic observables. Each source generates a 
characteristic pattern of relations between different observables, rooted in its field content and 
its transformation properties under chiral symmetry. The observation of such pattern in the 
current generation of, for example, nucleon and nuclear EDM experiments would then effectively 
pinpoint the dominant ff mechanism at the QCD scale [131 E31 EH ESI [361 E3) if one exists. 
Once this is known, the next step would be to trace the dominant effects up to the electroweak 
scale. 
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On the other hand, the formalism we develop can be easily adjusted to specific extensions 
of the SM. Any low-energy observable can be obtained by combining hadronic contributions 
from each separate source. Once the values of £0,3, <5o,3> w, £, and cr^g in a given model, and 
their running from Mt to Mqcd> are known, then the relative importance of the interactions 
constructed in the next sections can be reassessed to accommodate, for instance, a large hierarchy 
between these parameters. 



4 SO (4) chiral framework 

Below the QCD scale quarks and gluons are no longer convenient degrees of freedom to describe 
strong interactions, which can be rewritten in terms of hadronic EFTs. The implications of the 



Lagrangian in Eqs. (30), (34), (38), and (54) for the interactions among pions and nucleons 
(and delta isobars, since they are not much more massive than nucleons) at low momentum 
Q ~ m-n <C Mqcdj where m n ~ 140 MeV is the pion mass, are described by chiral EFT, an 
extension to arbitrary number of nucleons of xPT. At such momenta, pions must explicitly be 
accounted for in the theory, while other mesons can be integrated out. At lower momenta, 
Q -C m n , even pions (and delta isobars) can be integrated out, the corresponding EFT being 
called pionless EFT. This EFT finds applications in loosely bound nuclei, where it is much 
simpler to deploy than chiral EFT for it involves only short-range inter- nucleon interactions. 
However, the constraints from chiral symmetry, and consequently some predictive power, are 
lost. The form of the relevant pionless EFT interactions can be read from the following by 
discarding pion fields. 

The special role of the pion is a consequence of the approximate invariance of the QCD 
Lagrangian under the chiral symmetry SUl(2) x SUr(2) ~ SO (4). Because it is not manifest 
in the spectrum, which only exhibits approximate isospin symmetry, chiral symmetry must 
be spontaneously broken down to its diagonal, isospin subgroup, 5£/y(2) ~ 50(3). From 
Goldstone's theorem, one expects to find in the spectrum three massless Goldstone bosons 
that live on the "chiral circle" 5 3 ~ 50(4)/50(3). The explicit breaking of chiral symmetry 
gives pions a mass, but the way interactions at the quark/gluon level transform constrains 
pion interactions. In this section we show how these constraints arise in the CP-even sector 
for one particular choice of pion fields, and in the next section apply this technique to -fjT 
interactions. (Because observables are independent of the choice of fields, the method we use can 
be straightforwardly reproduced for other choices without affecting the physics.) Generalization 
to 5f7^(3) x SUr(3) is possible, but the SUl(2) x SUr(2) case is best adapted to nuclear physics, 
on which scale the strange quark is not particularly light. The method used in this section stems 
from Refs. [251 EH- 

We parametrize the chiral circle with stereographic coordinates [25J, whose dimensionless 
fields we denote by an isovector field £. We can identify these degrees of freedom with canonically 
normalized pion fields tv = F^C-, where F n ~ 186 MeV, the pion decay constant, is the diameter 
of the chiral circle. Such fields transform in a complicated way under chiral symmetry. However, 
a pion covariant derivative can be defined by 

Dfi"K = D~ 1 d^7v, (55) 

with 

D = l + ^ 2 , (56) 
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which transforms under chiral transformations as under an isospin transformation, but with a 
field-dependent parameter. Similarly, we can use an isospin-1/2 nucleon field N = (p n) T that 
transforms in an analogous way, and a nucleon covariant derivative, 

V^N = (dp + ■ 7T x D„tt) N. (57) 
We define P t through NV^ = VN and use the short-hand notation: 

nV^ = nV ± V^n, nV^V^ = nV^V u + V^V^n ± V^nV v ± P^TjP^. (58) 

Covariant derivatives of covariant derivatives can be constructed similarly, for example 



{V^D u ir)i = (dpSij - ^e ikj (iv x D^ir)^ D u 7Tj 



(59) 



For simplicity we omit the delta isobar here, but one can introduce an isospin-3/2 field for it 
along completely analogous lines joT) . 

Since nucleons are essentially nonrelativistic for Q <C ttin, the nucleon mass, we work in the 
heavy-baryon framework [31] where, instead of gamma matrices, it is the nucleon velocity 
and spin (S = (a/2,0) in the rest frame v = (0, 1)) that appear in interactions. Below we 
use a subscript _L to denote the component of a four-vector perpendicular to the velocity, for 
example 

V\ = V 1 - v»v ■ V. (60) 

The interactions we construct are manifestly invariant under rotations and translations, but not 
under Lorentz boosts. Nevertheless, Lorentz invariance imposes non-trivial constraints on the 
interactions in the effective Lagrangian and on their coefficients [42J. 

The first step in describing QCD at low energies is to construct the most general Lagrangian 
that transforms under the symmetries of QCD in the same way as the QCD Lagrangian itself. 
Along with this, one needs a power-counting scheme so that interactions can be ordered according 
to the expected size of their contributions. The Lagrangian contains an infinite number of terms 
that we group using an integer "chiral index" A [23] and the (even) number of fermion fields /, 

oo 

£ = EE4 A) - ( 61 ) 

A=0 f/2 

We will restrict ourselves here mostly to / < 2, leaving a discussion of the case / > 4 for Sec. 

oraj 

The technology for constructing the Lagrangian is well known, see, for example, Ref. [25]. 



When we neglect C e ( 34 ) , £4 ( 38 ) , and Cq ( 54 ) , the EFT Lagrangian includes all chiral-invariant 
interactions made out of D^tv, N, and their covariant derivatives. In this case, / < 2 interactions 
have chiral index [23] 

A = d + //2-2>0, (62) 

in terms of the number d of derivatives (and powers of the delta- nucleon mass difference). The 
coefficients of the effective operators, the so-called low-energy constants (LECs), cannot yet be 
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calculated directly from QCD, but they can be estimated using naive dimensional analysis (NDA) 
|55} I20|. in which case the index A tracks the number of inverse powers of Mqcd ~ 2irF 7T ~ 1.2 
GeV associated with an interaction. (Note that since NDA associates the LECs of chiral- 
invariant operators to g s /4:ir, for consistency one should take g s ~ 4-7T.) For the purposes of the 
present work, we need explicitly only the leading CP-even interactions, 

4°/<2 = -D^tz + nUvV- 2 -^S»t ■ N, (63) 

where qa — 1-267 is the pion-nucleon axial-vector coupling. At this order the nucleon is static; 
kinetic corrections have relative size 0(Q/Mqc~d) and appear in D- 1 '. 

The formalism to include chiral-symmetry-breaking operators in the SU(2) x SU(2) 
Lagrangian has been developed in Refs. |25LI54j . Operators that break the symmetry as compo- 
nents of chiral tensors can be obtained by rotating operators constructed with covariant fields 
^ such as nucleon fields, and nucleon and pion covariant derivatives, 

O^... w (7r,*) = E aa /(7r)i?^(7r)---^(7r)O a / (8 /...^(0,*), (64) 

where the chiral rotation R is given in stereographic coordinates by 

I °V D^f DFk \ 

R a( s(7v)=\ _ 2n . lfiTT^ • (65) 

V DF n D y Fj) J 

In the EFT these operators generate interactions, now involving 7r directly, that transform as 
tensors and their tensor products. The strengths of these interactions are proportional to powers 
of the symmetry-breaking parameters, other dimensional factors being estimated by NDA. 



The most important chiral-breaking term is the m term in Eq. (38). It generates in the 
effective Lagrangian all possible interactions with the structures S4 , S^S^, etc. with strengths 
proportional to, respectively, mr(6), [fhr{6)) 2 , etc. In a similar way, one can incorporate in the 



EFT [54, 56j the em term in Eq. (38), which leads to isospin violation as -P3 and its tensor 



products. The most important terms are, omitting a constant irrelevant for our purposes, 

r (o,2) 1 / 2 . Am^\ 2 5ml 2 

4,/=0 = ~2D { m * + * + 2~D^ (66) 

and 

C»__ 2 - A« (l - iViV + *%Lff („ - *S>r ■ .) N. ,67) 

Here the dominant contribution to the pion mass is m\ = 0(r(9)fhMQco)j its correction 
Am 2 = C(m^./MQ CD ), and the nucleon sigma term Amjy = \m\ /Mqcd), as the respective 
interactions have the structures of S4, S4 Cg) £4, and 54. In addition, the quark-mass contribu- 
tion to the nucleon mass splitting is <5m/v = 0(r(#) _1 em) = 0{r(ff)~ 2 em\j 'Mqcd) and to the 
squared pion mass splitting dm 2 = 0(5m 2 N ), as the respective interactions have the structures 
of P3 and P3 CS> P3. When we are interested in processes with typical momenta Q ~ m^, it 
is convenient to trade m by to^/Mqcd in all chiral- variant terms, which then have strengths 
proportional to powers of m\ times appropriate powers of Mqcd- Since e ~ 1/3 [57], we choose, 
for simplicity, to count it as 0(1). Since by all evidence 9 is small, we also take r{6) as of 0(1). 



15 



The power counting of / < 2 interactions, Eq. (62), is straightforwardly generalized by defining 
d to count powers of as well |23j . 

Other isospin-violating hadronic operators and interactions with the photon field come 
from C e , Eq. (34). For applications to EDMs, it is important to construct electromagnetic 



operators in which hadrons interact with soft photons (with momenta below -Mqcd) m a gauge- 
invariant way. We can minimally couple charged pions and nucleons to the photon by modifying 
their covariant derivatives, 



1 

D 



(d^ij + eAfj,e 3i j)irj, 



<9 M + r • (tt X D^ cm 7v) + —A^ (1 + T 3 ) 



N. 



(68) 



For brevity, in the following we omit the label "em". In addition, we can couple the photon 
through the field strength F^ u , in operators that transform as 7 M /6+T^/2 or its tensor products. 
The index A defined in Eq. (62) can be generalized to label electromagnetic operators, by 



enlarging the definition of d to count also the number of photon fields, which, having dimension 
one, require compensating powers of Mqcd in their coefficients. Finally, integrated-out hard 
photons give rise to purely hadronic operators that also transform as tensor products of + 
T34/2 [M], and are proportional to powers of (e/47r) 2 = a e m/4vr. The most important of these 
operators is the electromagnetic contribution 5m^ to the squared pion mass splitting, 



(i) 

em,/=0 



5ml 
2LP 



-1), 



(69) 



with 5m^ = 0(a era MQ CD / Ait). Since numerically a em /4vr ~ (to^/Mqcd) 3 , we assign d = 3 
for each power of a em /4ir [53]. This counting reflects both the smallness and the presumed 
electromagnetic origin of the pion mass splitting. 

The remaining dimension- four term, the to* sin term in Eq. (38), which is PP, can be 
treated in the same way. This was done in some detail in Ref. [32], where the implications 



of its P4 structure were discussed. Perhaps the most important point is that the hadronic PP 
interactions are intrinsically linked to isospin-breaking operators that transform like P3, so that 
the PP LECs can be inferred from isospin-breaking couplings, when the latter are known. 
In the rest of this paper we discuss the construction of the PP interactions stemming from Cq, 



Eq. ( 54 ) , and their power counting. The PP chiral Lagrangian can be constructed by writing 
down all terms that transform in the same way under Lorentz, P, P, and chiral symmetry as the 



terms in Eq. (54). Because such interactions are very small compared to CP-even interactions, 



we restrict ourselves to low-energy operators involving at most one power of the dimensionless 
coefficients <5o.3, <5o,3) w, £, and cr^g, even though it is straightforward to consider mixed effects. 
As discussed above, these coefficients are model dependent, and we construct the low-energy PP 
Lagrangian for each source separately. The form of the low-energy interactions is determined 
by the way a source breaks chiral symmetry. 

The gCEDM and the two four-quark operators in Eq. (54) with coefficients ImSi § do not 



break chiral symmetry and are therefore SO (A) scalars, which implies that they induce the same 
chiral Lagrangian and cannot be separated on the basis of low-energy experiments alone. From 
now on we refer to these as chiral-invariant sources (xISs) and use the symbol w (I) to denote 



collectively the dimensionless constants w and (Ti g (the invariants Iyy, Iqq , and J^) in Eq. (48) 
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(Eq. m): 



{w,ai,a 8 } 



w, 



I. 



(70) 
(71) 



Pion interactions stemming from xISs involve the pion covariant derivative, unless I is in a 
tensor product with a chiral- variant tensor (i.e., a product of Ss and Ps). 

In contrast to the xISs, the other sources (qEDM, qCEDM, FQLR) break chiral symmetry 
explicitly and can bring in pion non-derivative interactions by themselves, as is the case for 9. 
The specific dependence on n depends on the source. The qCEDM consists of two independent 
components, the isoscalar and isovector qCEDMs. They transform, respectively, as the fourth 



and third components of the 50(4) vectors V and W defined in Eq. (51). Since the QCD 9 



term transforms as the fourth component of an SO (4) vector as well, the 9 term and isoscalar 
qCEDM generate identical chiral operators (but with different strenghts). 

By themselves, the FQLR operators have the most complicated chiral structure, transforming 
as the 34-component of the symmetric tensors X^ and X^ in Eqs. (52) and (53). As for xKs, 



both operators induce the same chiral Lagrangian, but their entanglement is not a real issue 
since they depend on the same dimensionless parameter £. From now on we use the symbol X 
to denote X^ and X^: 



{x^,x^} -> X. 



(72) 



The structure of the qEDM, transforming as the fourth and third components of the 50(4) 



vectors V and W in Eq. (50), resembles that of the qCEDM, but the similarity is deceiving due 



to the photon. When the photon is integrated out to produce purely hadronic operators from 
qEDM, one needs to take the tensor product with the 1^/6 + from Eq. (34). Conversely, 



to produce operators with a soft-photon the qCEDM requires (as do FQLR and xISs) a tensor 
product with I M /G + T^/2. In both cases, the extra T^ A produces in general interactions of more 
complicated structure. 

Note that for each chiral-variant source, there are associated quark-gluon operators that 
do not violate P and T. For example, a qEDM is associated with a quark magnetic dipole 
moment, since they make the same 50(4) vector, Eq. (50). Analogous statements hold for 



qCEDM in Eq. (51) and for FQLR in Eqs. (52) and (53). The situation is similar to the 9 



term, which is related to the quark mass splitting [32J. However, for dimension-six sources these 
relations are less useful, since the corresponding OP-even interactions are buried among much 
larger, dimension- four interactions. 

The interactions stemming from the dimension-six fp sources can be organized according to 



a chiral index analogous to Eq. ( 62 ) , with the only difference that the coefficients of low-energy 



interactions must contain two inverse powers of the high-energy scale My, which replace two 
inverse powers of .Mqcd- The powers of Mqcd in a coefficient are therefore counted by 



A 6 = d + //2-4, 



(73) 



where d is the number of derivatives, and powers of the quark mass and e as described above. We 
will find that for the qEDM, qCEDM and xISs A 6 > -2, while for the FQLR operator A 6 > -4. 



The suppression by at least one power of a e 
hard photons essentially irrelevant. 



i/47r renders the purely hadronic operators from 
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The question now is to what order, for each separate source, do we need to construct the 
f JT chiral Lagrangian. This obviously depends on what observable one wants to calculate. The 
most likely candidates are -fjT electromagnetic moments and form factors of light nuclei, which 
require the use of the most important ff 1 pion, pion-nucleon, photon-nucleon, and inter-nucleon 
interactions. In order to compare our EFT approach to more traditional approaches, where 
PJT is implemented through three non-derivative pf pion-nucleon interactions, we construct the 
pion-nucleon Lagrangian for each source up to the order where all three pion-nucleon interac- 
tions appear. As found in Refs. |35} [36], for most sources the pf nucleon-nucleon (-photon) 
interactions are subleading with respect to ff one-pion exchange, so in the nucleon-nucleon 
sector we only construct LO operators. We construct the electromagnetic sector until, for each 
source, we find the momentum dependence of the nucleon EDFF, as the linear term in the 
square of the momentum transferred is a contribution the Schiff moment, which is important 
for the evaluation of atomic EDMs |28| . Since operators with two explicit photons give small 
contributions even to atomic EDMs [58] , we do not construct explicitly here operators with more 
than a single soft photon. 

We start in Sec. [5] with the construction of operators in the purely mesonic sector (that is, 
with / = 0). We then tackle the pion- and photon-nucleon sectors (/ = 2) in Sees. [6] and [8j 
respectively. We extend our analysis to include the nuclear sector in Sec. [TUj 



5 Pion sector 

In contrast to CP-even dynamics, p£ can generate interactions with an odd number of pions, 
starting with pion tadpoles that represent the disappearance of a neutral pion into the vacuum. 
Tadpoles cause the vacuum to become unstable because it can create neutral pions to lower its 
energy. In the 6 case, imposing the condition of vacuum alignment at the quark level in first 
order in the quark masses |44] ensures that the pf breaking is through the SO (A) vector P 



in Eq. (39) and that tadpoles are absent up to that order. For the dimension-six sources, the 
more complicated chiral structure also leads to tadpoles. Since the tadpoles are small, they can, 
in principle, be treated in perturbation theory, but a more convenient way of handling them is 
to rotate them away by performing a field redefinition [32] . After we consider the two leading 
orders for each source, we discuss the effects of tadpole extermination, which has important 
consequences for other sectors of the theory. 



5.1 qCEDM 

As it was found in Ref. |32] , without any nucleon fields it is not possible to construct an operator 
that transforms as the fourth component of an SO (4) vector, so the isoscalar qCEDM has no 
effect without another source of chiral breaking. The isovector qCEDM, however, breaks isospin 
symmetry as (I3W3 and generates a pion tadpole with chiral index Aq = — 2, 

^q,f=0 9 2D ' 1 ' 



where, by NDA, the LEC scales as 



A ,- 2 , = ^3"^™) . (75) 
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The next operators appear two orders higher in the chiral expansion. First, there is a ^3^3- 
type operator built from two pion covariant derivatives. Second, there are contributions from 
the combined effect of the quark mass and the qCEDM, producing pionic operators with the 
same chiral properties as the tensor products (I3W3 ® 77164 and do V4 ® fhePs. Together, these 
are 



(0) 

q,f=o 



^3 / n n » n 1 a (0) ^3 (, 27T 2 \ 
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where 



01=O[S3 



The + in the scaling of A~ should not be taken literally, but as an indication that the LEC 
gets contributions from two sources. 



(76) 



(77) 



5.2 FQLR 

In the pionic sector, X34 also leads to a tadpole operator, but of different form than qCEDM's: 
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with coefficient 
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Equation ( 78 ) generates different three-pion vertices than Eq. ( 74 ). As we will see, this difference 



is important, because the elimination of the pion tadpole does not completely cancel the operator 



in Eq. (78); instead, it leaves some three-pion coupling behind. 



Two orders down we find terms with two derivatives and terms transforming as ImHX34 
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with the scalings 
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5.3 qEDM 

The leading pion operators induced by the qEDM transform as ^3^3 ® el 11 /6 or (I0V4 eT^/2, 
and are simply given by 

-Al 1 )^, (82) 
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where the LEC scales as 

Af> = (^(*+*)=^E). (83) 

Again the + in the scaling of the LEC should not be taken literally. Because of the a e m/47r 
suppression, we do not bother to go to higher order. 



5.4 xISs 

In the mesonic sector it is not possible to write down tf^p chiral-invariant operators. Operators 
in this sector can be constructed by combining the effects of xISs with chiral-symmetry breaking 
from the quark mass difference, which renders them identical in form to those of the qCEDM. 
Indeed, at lowest order dy/I <8> emP^ gives 



A '(-2) _A(-2) i? T 7r 3 

with 



L wj=o ~ A ™ 2D~> ' ' 



A<-> = (» E ^2j. (85) 
At Aq = 0, we get operators with two covariant derivatives and operators transforming as 



(86) 



where the LECs scale as 




$i = 0[ we^ ] , A(°)=0( we^ } . (87) 



5.5 Tadpole extermination 

We have just seen that the transformation properties of the dimension-six sources cause pion 
tadpoles to appear in the ff mesonic Lagrangian. Since the coupling constant of the neutral 
pion to the vacuum is small compared to the pion mass, these tadpoles can be dealt with in 
perturbation theory, meaning that, for any given ^°JT observable at a given accuracy, only a finite 
number of neutral pions disappearing into the vacuum must be considered [32J. For applications, 
such as the calculation of hadronic [33] and nuclear |34^36j EDMs, it is however more convenient 
to eliminate the pion tadpoles from the mesonic Lagrangian, which can be achieved with field 
redefinitions of the form discussed in Ref. |32j . 
We define a new pion field Q' = tv' /F^ through 

where 

d! = 1 - C(l - C' 2 ) + 25^, (89) 
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and 



C = — (1 — cos ip), 



■ sin 99, 



(90) 



in terms of an angle (p. This transformation looks complicated, but has the nice property that 
the pion covariant derivative transforms simply as 



(91) 



with an orthogonal matrix 



(92) 



This simple transformation ensures that chiral- invariant operators, built from pion covariant 
derivatives, are invariant under the field redefinition in Eq. (88). This is not the case for 
operators that break chiral symmetry and involve the pion field directly. 

Before performing any rotation we summarize here the non-derivative, chiral-symmetry- 
breaking Lagrangian in the purely mesonic sector, including quark mass and ff operators 
in the first two orders (except for the very small qEDM operators, for which we keep only the 
least unimportant term). From Eqs. (66), ([7§, (re), @, @, (g4), and @, 
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We omit chiral-breaking operators generated by the electromagnetic interaction, which are not 
affected by the field redefinitions |32j. 

Our first goal is to remove for each source the dominant tadpole, that is, the terms with index 
A 6 = -4 for FQLR, A 6 = -2 for qCEDM and xISs, and A 6 = 1 for qEDM. By performing the 
field redefinitions in Eq. ( 88 ) with the small angle 
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(94) 



the dominant tadpoles are removed. From here on we keep only terms linear in ip. The effect of 
this field redefinition on the mesonic Lagrangian, apart from canceling the dominant tadpoles, 
is to modify the coefficients of the tensors: 
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in terms of the shifted LECs, 



A[ _ 2) , = _4?%^A<->= (W)^^), (%) 



mt 




ASP- = KV-^^^ = oUl + e>)<y (98) 

(Here and in the following we neglect terms of higher order than we need in the rest of the 
paper.) We can thus rotate away the leading tadpoles without introducing new interactions in 
the meson Lagrangian. The net effect of the rotation is only to change the dependence of the 
coefficients on the parameters 5q, 63, and e. 



The remaining tadpoles in Eq. (95) can be eliminated by a second rotation, now with an 
even smaller angle 

«> = (4e? + AS + A< R 2 " + Af + A<°><) , (99) 

which leaves us with 



tadpole " -w{ 1+ ri*D) 77 + 2^ 3 

"fS ( A - + A -' + ^ + 4 0) ' + ^ 0)/ ) • (100) 

Although the tadpoles are removed, residual f^P interactions involving an odd number of pions 
are left behind. In case of the qCEDM and xISs these terms carry a high chiral index and such 
multi-pion vertices contribute to observables at high order. In case of the FQLR, due to its 
complicated SO {A) properties, a three-pion interaction remains with a relatively low index. The 
consequences of this interaction are discussed in Sees. 7_A, 9j and 10. 5[ 



So far we did not discuss the operators with LECs These operators are not affected by 
the above rotations up to the order we are working. They induce changes in the pion kinetic 
term, which nevertheless can be removed by another pion-field redefinition. 



6 Pion-nucleon sector 

For most applications, we can ignore effects in the pion sector, once the tadpoles have been 
removed. Tadpole extermination does leave traces, however, in operators involving nucleons. In 
this section we construct the ffi pion-nucleon operators, incorporating the new terms arising 
from the field redefinitions needed to remove the tadpoles, pion-nucleon interactions are 
very important for nucleon and nuclear EDMs. 
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6.1 Tadpoles strike back 



It is extremely convenient to supplement the pion field redefinition (88) with a nucleon field 
redefinition |32| 

N = U'N', (101) 

in terms of the unitary matrix 



'd! 



U' 



(102) 



where d', C, and S are given by Eqs. (89) and (90). This redefinition ensures simple transfer 



mation laws for covariant objects such as nucleon covariant derivatives and nucleon bilinears, 



V.N = U'V'N', 



NN = N'N', 



NttN 



O'ijN'rjN', 



(103) 



where O is the orthogonal matrix in Eq. (92). The properties (103) extend to the nucleon sector 



the invariance under field redefinition of chiral-invariant operators, built from nucleon fields and 
covariant derivatives. 

The two rotations that allow elimination of the pion tadpoles in leading orders affect the 
chiral-breaking terms involving nucleons. Most importantly, the nucleon sigma term and mass 



splitting in Eq. (67) become, after the rotation of Eq. (94) 
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The terms linear in (p contribute to the -fjT pion- nucleon interactions, in addition to those 
generated directly from the chiral structure of the PT dimension-six sources. Higher-order pion- 



nucleon operators come from rotation with angle ip' (99) and/or higher-order chiral-breaking 
terms in the CP-even Lagrangian. 

We refrain from giving the boring details here. Instead, we turn to the construction of the 
PJT pion-nucleon interactions from the chiral structure of the various sources. We indicate below 
the changes that arise from tadpole extermination. 



6.2 qCEDM 

Since the qCEDM breaks chiral symmetry, it gives rise to non-derivative pion-nucleon interac- 
tions. Already at chiral order Aq = —1, the structures doV^ and d^W^ give 
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Comparing Eqs. (|105|) and (|104[) shows that tadpole removal leads to shifts 
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While g\ gets a change of 0(1) that nevertheless does not affect its scaling, go now depends on 
the isovector qCEDM as well. In the following we absorb the terms oc A^ in g^. 

Equation ( 105 ) shows the first important difference between from dimension-six operators 
and the QCD 6 term, namely, the presence of the ffi isospin-breaking interaction with LEC g\ 
at leading order in the / = 2 Lagrangian. This fact is particularly relevant for the f^T moments 
of the deuteron [591 EU [36] . The 9 term also generates this interaction but it is suppressed by 
(ttItt/Mqcd) 2 compared to leading-order interactions [32J. 

Increasing Aq by one, we find operators with the same chiral structures but one covariant 
derivative, 

4°/= 2 = ^(K-V^NSm + ^nsV^.NS'TN 

+ k (? ?3 ~ yB) mrxv ' D7T)tN ' (109) 



where the LECs are 
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Increasing A6 by another unit, the chiral structures proliferate significantly, and so does the 
number of different interactions. With still the same chiral structure as before, but two covariant 
derivatives, we find 
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The scaling of the LECs Q and £j is given by 



2 \ II 

ra t \ - ,„ / ~ mi 



Ci = O[S n , & = O U 3 „ * . (112) 

Sl ^ M|,Mqcd y ? [ 3 M^M QCD J { ' 

Reparametrization invariance relates some of these operators to operators of lower order, 

7 _ 7 90 7 _ 9A9o h 

U — C,6 — ~ — o", <s8 — o ' 

2m^ m^f ?n jv 



91 7 03 7 9A9\ 

In these subleading pion-nucleon Lagrangians, the effects of the field redefinition on operators 
that transform as W3, like those with LECs ^2 and 03 in Eq. (109) and in Eq. (Ill), can 
be absorbed in a redefinition of the coefficients, whose scaling is still given by Eqs. ( |110 ) and 
(112). For operators that transform as V4, the scaling of the coefficients is modified, and they 
get a contribution from the isospin-symmetry-breaking Lagrangian at order A = 1, 2 (listed, for 
example, in Ref. |32j). Schematically, 

A ^S = o ((«.+*) 5 i^)- a") 

Still at order Aq = 1, interactions arise from the combined effect of the qCEDM and the 
QCD mass terms. The resulting operators transform as (mS^ — emP^) <g) (— d^V^ + ^3^3). They 
generate 

rW 5 90 A 2tt 2 \ - 5g! ( 2tt 2 \ - 

£ ^ 2 = ~i^V-^) Nt ^ n -f^{ 1 -^)^ nn 

112 -ir 3 N (t 3 - J^-r • tt^) AT. (115) 



The first two terms are corrections to the isoscalar and isovector couplings in Eq. (105), 
from which they differ only by terms with three or more pions. The corrections scale as 
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For most practical purposes 5 go and 5g\ can be absorbed into go and g\. The third operator 
is the most interesting, as it is the first contribution of the qCEDM to the isospin-breaking 
pion-nucleon interaction tt^NtsN. It has the scaling 

^°( eSs ^)' (U7) 

Thus, two orders above lowest all three possible non-derivative pion-nucleon f*p interactions 
receive a contribution from the isovector qCEDM. Just like the 9 term |32j, the isoscalar qCEDM 
generates 52 as well, but this requires a photon exchange so that 32 is suppressed by a eTa /4ir and 
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enters three orders above lowest. This provides the second difference between isovector qCEDM 
and term (and isoscalar qCEDM). However, since for all sources 52 enters in the subleading 
Lagrangian, this difference is of little phenomenological interest. 

Tadpole elimination affects subleading non-derivative couplings in two ways. First, the 
elimination of subleading tadpoles shifts the coefficients go and g~\ in Eqs. (107) and (108), 
5o,i -> 50,1 + tf'go.l, witri 
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(118) 



contributing to the Aq = 1 Lagrangian. Second, the operators in Eq. (115) receive corrections 



from the transformation of CP-even, chiral-breaking operators that contribute to the chiral 
Lagrangian at order A = 3. These operators are proportional to two powers of the quark 
masses, and transform as the tensor products {mS^ — rheP^) (g) (77164 — meP-s). The rotation 
needed to eliminate the leading tapdoles causes them to generate corrections to the coefficients 
5 go 1 and 52- Now 5 go also depends on eS^ and 5g\ receives a £ 2 5s correction, 



^5o -> 5g = 0\ [5o + eS, 
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while the scaling of 52 is unchanged. 
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6.3 FQLR 

The FQLR operator also can generate non-derivative pion-nucleon couplings. Pion-nucleon 
operators start at Aq = —3, and at this order there is only one operator with the structure 
of X34, the isovector pion-nucleon interaction. It appears as if the FQLR generates only 5i at 
LO, in stark contrast with the 9 term and qCEDM which generate, respectively, go and 5o,i & t 



LO. However, the removal of the tadpole in Eq. (78) effectively causes the appearance of the go 
interaction at Aq = —3. At this order, 
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where 



arises directly, and 
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stem from tadpole removal via Eq. (104). One can eliminate Aljj 4 ' 1 and write 
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A lattice evaluation found Stun = 2.26 MeV [60j, while values for Amjy range between 45 and 
60 MeV [61] . Therefore, for FQLR go is actually only a few percent of g[. 

At higher orders FQLR generates some interactions that resemble those of the isovector 
qCEDM, although with a more complicated structure, plus some new interactions. At Aq = — 2 
and Ag = — 1 we find interactions similar to those with LECs ^23 in Eq. (109) and £j in Eq. 
(Ill), but with the replacements 
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The additional interactions are 
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and 
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^2,3,4 = O 



6 = o e : 



M, 



QCD 



7 1 



The RPI relations are the same as in Eq. (j 1 1 3[) , with the additional relation 
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Just like for qCEDM, tadpole removal induces the subleading isospin-conserving f^p operators 



with LECs ft\ in Eq. (109) and Q in Eq. (Ill), but with scalings 
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Again, two orders above the lowest order, operators appear due to an insertion of the quark 
mass (difference), 
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where the LECs scale as 
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Although the operators in Eq. (130) appear to be very complicated, if we ignore terms with 



three or more pions the Si 291 and Sgo couplings are contributions to the standard non-derivative 
pion-nucleon interactions. Just as for the isovector qCEDM, the third pion-nucleon coupling g<i 
comes in two orders higher than 50,1 • 

As in the case of the qCEDM, the elimination of the tadpoles modifies Eq. (130). The 



- c_2) - (—2) _ _ 

elimination of the subleading tadpoles Al R1 and A^ R2 shifts Sgo and 5\gi, but does not modify 
their dependence on £ and e. On the other hand, the effect of the elimination of the leading 
tadpole on chiral-breaking operators proportional to two powers of the quark masses m and fh e 



is to generate a contribution identical to Eq. (115), with the coefficients replaced by 
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6.4 qEDM 

The pion-nucleon interactions originating from the qEDM arise from the tensor product (— doV^+ 
(I3W3) ® e(/ M /6 + and thus have a rich chiral structure: 
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The qEDM also generates tadpoles, as given in Eq. (82), which can be removed in the same 



way as for the other sources. This removal alters the e dependence of go , in a way that can 
be summarized by replacing (So + S3) with (So + ^3) (1 + e) in Eq. (134). Since, as Eq. (134) 
shows, 50,1 already receive contributions from both the isoscalar and isovector qEDMs, the added 
contributions are not particularly relevant. 

Note that all possible non-derivative pion-nucleon interactions appear at the same order. 
(This only holds for the isovector qEDM, since for the isoscalar qEDM g~i is suppressed.) Only 
the qEDM has this property: for all other ff dimension-six sources 52 is suppressed with respect 
to go and/or <?i- However, the a em /47r suppression makes all the qEDM pion-nucleon couplings 
pretty unimportant. 



28 



6.5 xISs 

In the pion-nucleon sector, no chiral-invariant fp operator can be constructed with zero or one 
covariant derivative. The first operators therefore start at Aq = — 1 and have two covariant 
derivatives or one insertion of the quark mass (difference) : 

</l 2 = ^(V^ ± D^).NtN+^(D,7vxv.Dtt).NS^tN 



with scalings 
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It is interesting to notice that rotation invariance alone would allow a two-derivative operator of 



the form ^(D v tt) ■ Nt[S^ , S U ]T>^N in Eq. (135). However, this operator is not RPI by itself, 



and its variation cannot be absorbed by any other operator in the leading-order Lagrangian. 
RPI, therefore, forces £4 to vanish in leading order. The operator (4 appears at NNLO, and it 



is linked to go by RPI, as in Eq. (113). 



Here the tadpole rotation induces the shifts 
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9 i -> g x + Amjv— = O I ew— ^ J , (138) 

which are less interesting, because the dependence on w is not fundamentally changed. As 
before, we absorb contributions from tadpole elimination in gi 



From Eq. (135) we see that the xISs, just as the qCEDM and in contrast to the 9 term, 
induce g\ at the same order as go. Differently from all the other sources, the xISs also generate 
two- derivative operators of the same importance as go and g%. 

One order higher we find various types of interactions: chiral-invariant operators with three 
covariant derivatives; chiral-symmetry- breaking operators with one power of quark masses and 
one covariant derivative; and electromagnetic operators coming from the tensor product d\yl ® 
e(J^/6 + T^/2) ® e(/ M /6 + T^/2). The chiral-invariant operators are 

•(0) _ **l<n /-... rwt. at^ at , ix 2 
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where 
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while RPI gives 



The quark-mass insertions give rise to 
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The removal of tadpoles does not affect the chiral-invariant operators in Eq. (139), while it 
slightly modifies the e dependence of f3\ in Eq. (142), in a way analogous to Eq. (137). 



Finally, integrating out a hard photon gives the single operator 
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We see that at this order the first contribution to the non-derivative pion-nucleon interaction g 2 
appears. It scales as 

lem ^QCD \ 
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With the usual counting a em /47r ~ m^/Mq CD , this interaction is suppressed by one order 
compared to the 50,1 interactions. 



7 The pion-nucleon form factor 

An important element in the evaluation of hadronic and nuclear EDMs is the pion-nucleon 
coupling. In this section, we summarize the pion-nucleon interactions for the different 
sources we have considered, by calculating the ff PNFF with the Lagrangian derived above. 
For convenience, we consider the Lagrangian after tadpole extermination. 

We consider the three-point Green's function for an incoming (outgoing) nucleon of momen- 
tum p^ (p' '^) and an outgoing pion of momentum = p^ — p' ^ and isospin a. We take the 
incoming and outgoing nucleon to be nonrelativistic and on-shell, so 

P = 7, Am N =F -TT- + • • • , P = 7, Am N =F -^r~ + ■■■, (146) 

2m n 2 Zuin 2 

where the — (+) signs holds for protons (neutrons) and the arrow denotes vectors in three- 
dimensional Euclidean space, p^ = (p°, p), p^ = (p°, —p). It is clear that p° enters at one order 
higher than p. The Green function for on-shell nucleons can be parameterized by three form 
factors, corresponding to three different isospin structures, 

V a (q, K)=~ [F^q, K)r a + F 2 (q, K)5 a3 + F 3 (q, K)5 a3 r 3 ] , (147) 
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in terms of the functions F± 2 3 of and = (p^ +p ,fl )/2. In what follows, we give F\, F 2 , 
and F3 for each of the dimension-six PT sources. We calculate the PNFFs up to the order where 
each of them gets a non-vanishing contribution. As derived in Sec. [6j this implies going to 
NNLO for qCEDM and FQLR, and to NLO for xISs. For qEDM we need the LO contributions 
only. 

At LO and NLO, for all sources but FQLR contributions arise exclusively from tree diagrams. 
At NLO for the FQLR and at NNLO for the other sources loops appear. We use dimensional 
regularization in d spacetime dimensions, which introduces the renormalized scale \i and 

L= -^- d - 7 E + ln4vr, (148) 

where je — 0.557 is the Euler-Mascheroni constant. At this order we use q° = 0. 



7.1 qEDM 

The PNFFs from the qEDM are very simple to the order we are interested in, since all three 
PNFFs appear at the same, leading order: Ag = 2. We read off from the Lagrangian directly, 

Fi=g , F 2 = gx, F 3 = g 2 . (149) 

Since these interactions contain suppression factors of a era /4ir ~ 7t4/Mq CD , they are most 
likely not important in calculations of nuclear EDMs, which are dominated by short-range 
contributions to the nucleon EDM |36j . 



7.2 x ISs 

For the gCEDM and chiral-invariant four-quark operators, the PNFFs receive contributions at 



LO and NLO from interactions with LECs, respectively, go, gi, and £i in Eq. (135), and g 2 in 



Eq. (144). The contributions to the PNFFs can be read off easily, 



Fi(q) = go-g2-(if, (150) 
F 2 = gi, (151) 
F3 = 92. (152) 

We conclude that, for the chiral-invariant PT sources, F± and F 2 appear at the same order 
and -F3 appears one order down in the Q/Mqcd expansion. Apart from contributing to F%, 
g 2 also contributes to F\. This is of no phenomenological interest, since F\ receives a larger 
contribution from g~Q. At LO F\ depends on the pion momentum because of the presence of 
the chiral-invariant operator C,\ . The consequences for the PT NN potential have been worked 
out in Ref. [35], where it was shown that (j generates a long-range PT potential, which can be 
accounted for by a redefinition of go, and a short-range PT potential, which can be absorbed in 
a nucleon-nucleon contact interaction. 



7.3 qCEDM 

In LO (Ag = —1) and NLO (Aq = 0) the PNFF arises from tree diagrams where the pion- 
nucleon vertices are, respectively, the two non-derivative interactions with LECs go and g\ in 
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Figure 1: One-loop contributions of relative order 0(m 2 r / \2ttF 1t ) 2 ) to the pion-nucleon form 
factors Fi(q,K) and F2(q,K). A nucleon (pion) is represented by a solid (dashed) line; the ff 
vertices from Eq. ( 105| ) are indicated by a square. The other vertices represent LO interactions 
from Eq. ( 63 ) . For simplicity only one possible ordering is shown. 



Eq. ( [105] ), and the single-derivative interaction with LEC (3 3 in Eq. ( [T09J ). At NNLO (A 6 = 1) 
there are further tree-level contributions from the two-derivative interactions with LECs £i, (4, 
Ce, £i) £.4, £5 an d £7 in Eq. (Ill), and the non-derivative interactions with LECs 5go, 5gi and g% 



in Eq. (115). Use of the RPI relations (113) shows that the PNFFs can be expressed in terms 
of go, §1, /S3, Ci) £i) $9o, Sgi and g~2- The operators with LECs go, Ci an d $9o contribute to F±, 
those with LECs gi, £1, Sgi to F2, and that with LEC 52 to F3. The operator with LEC /3 3 and 
its recoil correction contribute to the PNFF 



V a (q,K) 



fa 
F v 



K 



m N 



(153) 



When we use the nucleon on-shell conditions (146), this contributes to F\ and -F3. 



At NNLO there are also one-loop diagrams where the vertex is one of the two non- 
derivative interactions with LECs go and g\ in Eq. ( 105| ), other vertices coming from the LO 
chiral Lagrangian, Eq. (63). These loops are shown in Fig. [I] The structure of the diagrams is 
such that the momentum of the external pion never flows into the loop and the only scale in the 
integral is the pion mass. As a consequence, the diagrams do not yield any non-trivial momentum 
dependence and simply renormalize the LECs 5g' 1 and 5'go,i- We define the renormalized LECs 



$90 



Jo + S' 'g + 



$9i = $9i + S'gi + 



go rn- 

4 (27^)2 
91 ml 



mz 



(1 + 3g 2 A ) [L + 1 - log -£ - 2g\ 



5 - 9 5 i) ( L + 1 - log ^ ) + (Sg\ 



4 (2^)2 

With this definition the PNFFs for on-shell nucleons read, up to NNLO, 

1 



F\{q, K) = go {I 



F 2 (q,K) = g^l 



2m 2 N 



K A + S-[K xq 



+ 5g + Sm N (3 3 - Ciq 2 , 



1 



2m 2 N 



K 2 + S-[Kxq) ^Sgx-Ciq 2 , 



g 2 - 8m N /3 3 . 



(154) 



(155) 

(156) 
(157) 



We see that F\ and F2 receive contributions at the same order. Two orders down we find 
momentum dependence of these PNFFs, as well as the first static contribution to F 3 . The 
implication of the momentum dependence to the two-nucleon potential was discussed in Ref. 

EH. 
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Figure 2: One- loop contributions from purely mesonic f^f interactions to the pion-nucleon form 
factors. The square denotes the tf^f vertex from Eq. (100). The other notation is as in Fig. [T] 



7.4 FQLR 

The main part of the PNFFs from FQLR is very similar to those from qCEDM. After replacing 
the scalings of the LECs, one sees that contributions to F\ and F 2 start at LO (Aq = —3), and 
contributions to F3 at NNLO (Ag = —1); at this order Fi and F2 obtain analytic momentum 
dependence. 

Again, the diagrams in Fig. [TJonly renormalize LECs. The different chiral structure of the 
isovector vertices g\ and g[ in Eq. (120) only affects the first diagram in Fig. [TJ with the 
minor consequence of modifying the counterterm 5g\ with respect to Eq. ( 154[ ). The results 



in Eqs. (155), (156), and (157) give the dominant contributions to the three PNFFs, with the 



replacement g\ - 
$91 



gi+g v and 

5 

= £151 + $291 + 5'gi + - A (351 + g'x 



ml 



+ 



3<& 



+ J 



rn 



(2irF n 



L + 1 - log 



mt 



{2ttF 7T ) 



3 [L + 1 - log 



mt 



+ 2 



(158) 



However, this is not the whole story. As shown in Sec. 5.5, the elimination of the tadpoles in 



case of the FQLR leaves the three-pion vertex with LEC A^ 4 ^ in Eq. ( 100 ) with a lower chiral 



index than the dominant pion-nucleon interactions. The one-loop diagrams in Fig. [2] contribute 
to the PNFFs at NLO (Ag = —2) and add to F2 a non-analytic momentum dependence not 
present for the qCEDM. This PNFF becomes, instead of Eq. 



(156), 



F 2 (q,K) 



(5l+5i) 1 



2m% . 



K + S ■ (K x q 



+ 5gi 



m A ( " 4) 



2m„ 



where 



/ 2 (x) = l + 



Using Eq. (122) we can eliminate Ajj^ in favor of g[ or go . Precise measurements on the 
deuteron ff electromagnetic form factors, which depend strongly on F2 [59. 1621 IM] . could, in 



l + 2x 2 
2x 



arctanx. 



(159) 



(160) 



principle, measure this momentum dependence and separate the FQLR from the qCEDM. 



8 Electromagnetic interactions 

The experimental interest in EDMs bring ff electromagnetic interactions to the forefront. 
Some interactions of hadrons with soft photons can be obtained using the {/(l)-gauge covariant 
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derivatives (68) in existing operators. More interesting are the interactions that arise through 
the field strength F^, which we describe here. 

Since the pion has spin 0, we cannot construct an EDM operator in the / = sector. In 
contrast, there are plenty oiffi interactions in the / = 2 sector. 



8.1 qCEDM 

In the case of the qCEDM, operators containing the electromagnetic field strength have the chiral 
properties of the tensor product of the source and the CP-even electromagnetic interactions. 
In lowest orders we need only (-J0V4 + d^Wz) ® e(I^/6 + 3f 4 /2). 

Interactions start at Aq = 1 and transform as the fourth component of a vector, or as the 
product of a vector and an anti-symmetric tensor, 



'(1) 

"q,/=2 em 



-2N 



dn 1 



2^3 
FID 



T ■ 7T 



S^Nv u F, 



fJbV 



-2d[ ( 1 



2tt 2 

FID 



N 
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7T T3 — 7T37T ■ T I 



+ 



1 



F V D 

+C2 



— ^ ua/3 v a N 



F*D 

V a N (C T • 7T + Cl7T 3 ) SpN F^ u 

2 , o 



S^Nv u F, 



F^D 



T3 



FID 



I 7T T3 — 7T37T • T I 



SpNF, 



/J.V, 



(161) 



where the electromagnetic LECs scale as 



do,i> c o,i = O e (S + 5; 



mz 



mz 



M^Mqcd J ' 



C-2 




mz 



(162) 



The first term in Eq. (161) is a short-range contribution to the isoscalar nucleon EDM. The 



second and third terms both contribute to the isovector nucleon EDM, and they differ only by 
interactions involving two or more pions, such that separating them is practically impossible. 



The last three terms in Eq. (161) are pion-nucleon-photon interactions which play a role in the 



calculation of the deuteron MQM [331137] . 

As always, the isoscalar qCEDM generates the same interactions as the 9 term [32]. But, while 
in the pion and pion-nucleon sectors the isoscalar qCEDM generated only isoscalar interactions 
in LO, in the electromagnetic sector both isoscalar and isovector interactions appear, due to 
breaking of chiral symmetry by the quark electric charge. The main extra feature of the isovector 
qCEDM is the appearance of the c~2 interaction, which has a T37T3 structure. However, it is 
unlikely that this operator has any important physical consequences. 

The elimination of the leading tadpole modifies to coefficients in Eq. 
can be schematically summarized by the replacement 5q 



5 + e5 3 in Eq. (162) 



(161) in a way that 
Since most 



operators already receive contributions from both the isoscalar and isovector qCEDM these 
shifts are not particularly interesting. 

The Schiff moment of the nucleon [431 155] and EDMs of light nuclei |36j are dominated by 
f^p pion-nucleon interactions. The LO nucleon-photon operators have little impact on these 
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observables, so we do not construct operators with higher chiral index, which would be even 
more suppressed. 



8.2 FQLR 

The interactions stemming from FQLR are similar to those from the qCEDM. Already at LO 
(Ag = —1) FQLR generates operators with a more complicated pion structure due to the 
X34 (g> tensor product in ImS X34 (g> e(/ M /6 + T^/2). However, apart from terms with two 

or more pions, the Lagrangian is identical to Eq. ( 161| ): 



^LR,/=2 em — ^q,f=2 em + ' • • ' (163) 

but with the scaling 

d ,i, d[, co,!, c 2 = • (164) 

Tadpole extermination brings in no important new features. Like for the qCEDM, operators 
with additional pions and/or higher chiral index have no obvious phenomenological interest. 

8.3 qEDM 

Since the qEDM contains a photon field, it yields nucleon-photon operators that transform like 
itself, namely as the third and fourth components of, respectively, the vectors V and W in Eq. 



(50). At LO this generates 

2 



C {1 ] , = -2N 

q,f=2 em 



- / 2tt 2 \ - / 2tt 3 \ 



S»Nv v F^ 



+— e^v a N(c r ■ 7v + CLn 3 )Sf}N Fur, (165) 



with 



2 \ 1 

vnt \ t „ / „ rat 



In contrast to the qCEDM, the isoscalar (isovector) qEDM generates isoscalar (isovector) nucleon- 
photon interactions, since the symmetry properties of the qEDM need not be mixed with chiral- 
symmetry breaking due to the quark charge. 

In the case of the qEDM, long-range physics propagated by pions is suppressed by powers of 
a em , and observables are dominated by short-range nucleon-photon interactions. Since the 



operators in Eq. ( 165 ) contribute to the nucleoli EDM only and not to the momentum-dependent 
part of the corresponding form factor, for the latter we need to construct electromagnetic oper- 
ators with higher chiral index. It turns out that momentum dependence arises only at NNLO, 
so that we need to construct the Lagrangians with Aq = 2,3. 

At Ag = 2, the ff electromagnetic operators that are not constrained by Lorentz invariance 
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contain at least one pion, 



(2) 

q,f=2 em 



-N 



m N 



27T 2 \ 



2vr 3 

+ d\ r 3 — t^tt • 7r 



2F- K rriND 



FID 



S^VK _NF, 



1 



J- 7T 



Pi 1 



2tt 2 \ 



2tt 3 



7T • D^7V 



Nv u F, 



+ ^-iV (p 2 7T • r + p 2 Tr 3 )Nv u d> i F tlv 
+ J^e°^V a N [p 3 (7V ■ Dpiv) + p 3 7T3 (^w) • t] iV F f 



III' 



+ ±_ e «^» R [^ Da7V + p 5WaU . D7V ) x t] 4 S^iV ^ e - 2 ^^j 



F 



llV-i 



where pi (pi) originate from the isoscalar (isovector) qEDM, with 



pi = ed 



mt 



mt 



(167) 



(168) 



In Eq. (167) we already incorporated RPI, and the first two sets of interactions are recoil 
corrections to Eq. (165). 

At Ag = 3, we encounter operators with two covariant derivatives or one insertion of the 
quark mass. Since it is unlikely that, at this order in the chiral series, operators containing 
pions are of any phenomenological use, we focus on terms without pions: 



C 
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q,f=2 em 



-2N (Sd + 5dm) S»Nv v F^ 
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(169) 



The first two terms are corrections to the isoscalar and isovector nucleon EDM due the quark 
mass, the third and fourth terms are relativistic corrections, and the operators with LECs S' and 
S[ are the first qEDM contributions to, respectively, the isoscalar and isovector Schiff moments. 
The dots denote multi-pion components of the listed operators, and other operators that start 
at one pion, which we neglect. The LECs scale as 



Sd = O\e (S + e5; 



Sd 1 = 0\e (s5 + 5; 



S' n = 0\ eSr 



mz 



S[ = 0[ e5 3 
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(170) 



Here tadpoles are much smaller than the nucleon-photon interactions, and their removal does 
not affect the operators above. 
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8.4 xISs 

The chiral-invariant sources generate electromagnetic interactions that transform as wl 1 
e(J"/6 + T 3 M 4 /2) with chiral index A 6 = -1, 



where the LECs scale as 



T3 + 



FID 



{ 7T37T • T — 7T T 3 I 



S»Nv v F„ v , 



(171) 



(172) 



Differently from qCEDM and FQLR, for chiral-invariant -fjT sources the short-distance EDM 
operators have the same chiral index as the leading pion-nucleon coupling, Eq. (135). Since 
the latter contributes to the nucleon EDM only via loops, it follows that, for xISs, the nucleon 
EDM is mainly determined by short-distance physics. 

A second consequence of the enhancement of short-distance vs. long-distance physics for 
XISs is that the nucleon EDFF does not depend on the momentum transfer at leading order in 
^PT. Momentum dependence only arises at NNLO, for which accuracy we need to consider the 
power-suppressed Aq = 0, 1 Lagrangian. We construct the complete chiral Ag = Lagrangian, 
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At Ag = 1 we only construct the operators that start without pions, since operators with 
pions are of little phenomenological use. We have 
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(175) 
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As for qEDM, the first operators (with LECs <5do l an d Sd^) are corrections to the isoscalar and 
isovector nucleon EDMs proportional to the quark masses. The next interactions are relativistic 
corrections to Eq. (171). Finally, the last operators (with LECs Sqi) are short-distance con- 



tributions to the first derivative of the nucleon EDFF, the Schiff moment. The scaling of the 
LECs is 



5d ,i = e(l + e)w 



5d[ = O 



ew- 



M 2 M qC D / 



Sn 



0,1 



O 



ew 
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(176) 

The leading nucleon-photon interactions in Eq. (171) are chiral invariant or transform as 



the 34 component of an antisymmetric tensor. In both cases, they are not affected by the field 
redefinition in Eq. (91). Only the NNLO operators <5do,i an d Sd'i, which are proportional to the 



quark masses, are affected by the elimination of the tadpoles, in a way that leads to the shift 
w — > w(l + e 2 ). 



9 Nucleon EDM from the FQLR 

The electromagnetic interactions of the previous section allow the calculation of EDFFs. In 
fact, the chiral Lagrangians from the qEDM, qCEDM, and xISs derived in this article have 
already been used to calculate the form factors of the nucleon and several light nuclei in 
Refs. \43> \ 133 1 134 ^ 136 ] 137] . However, in these references the contributions from the FQLR operator 
were omitted. In this section we update our previous work by giving the FQLR contributions 
to the nucleon EDM to NLO. Nuclear issues are discussed in the next section. 
The nucleon EDFF, following Refs. [291113], is decomposed as 

Fi{Q 2 ) = di- Sffi + Hi(Q 2 ), (177) 

where d% is the isospin-i component of the EDM, S\ is the corresponding Schiff moment [28], 
and Hi(Q 2 ) accounts for the remaining dependence on Q 2 = —q 2 > 0, q being the outgoing 
momentum of the photon. The EDFF of the proton (neutron) is Fq + F\ (Fq — iq). 

The LO pion-nucleon and nucleon-photon interactions from the FQLR are very similar 
to those from an isovector qCEDM, apart from interactions involving multiple pions, which 
contribute at higher orders. Therefore, the LO nucleon EDFF from the FQLR is of identical form 
as that from the isovector qCEDM calculated in Ref. [43], but, of course, with different scalings 
for the LECs. The LO nucleon EDFF gets contributions from the short-range interactions in 



Eq. (163) and from one-loop diagrams involving the LO isoscalar pion-nucleon vertex go m 
Eq. (120) (c/. Ref. [3] [28. 29J). At NLO, there are additional one-loop contributions, involving 
both go and the isovector vertex g\, and subleading couplings in the PT-even Lagrangian (c/. 
Refs. [311 HOI |33] ) . One has in principle to consider also the effects of the three-pion coupling 



A^pj 4 ^ in Eq. (100), which contributes to the nucleon EDM via two-loop diagrams. However, it 
is easy to see that all the two-loop diagrams vanish because of their isospin structure, so that the 
nucleon EDM at NLO stemming from the FQLR is identical to that generated by the qCEDM. 
In the evaluation of loop diagrams we use dimensional regularization, as in Sec. [7| 

The leading loop diagrams do not generate an isoscalar EDFF, which is therefore purely 
tree-level and static at LO. The first non-analytic contribution to the isoscalar EDM arises at 
NLO, and it has both go and gi pieces. At the same order one also finds the first momentum 
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dependence of the isoscalar EDFF, which is proportional to the hadronic part of the nucleon 



mass splitting, 5m^ in Eq. (67), and the isoscalar coupling go. The isoscalar EDFF to NLO is 
given by 
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such that /ig 1 ^ (x 



x 2 + 0(x 3 ) for x < 1. 

In contrast, for the isovector EDM loop diagrams do renormalize short-distance operators, 
in addition to generating a non-trivial momentum dependence. Instead of the nucleon mass 
splitting from the quark masses, there is a contribution from the electromagnetic pion mass 



splitting 5ml i R Eq. (69). To NLO, the isovector EDM is found to be 
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are defined so that they satisfy /i^(x <C 1) = x 2 + C(x 3 ). 

The dependence on the arbitrary scale /i from the loops is compensated by the counterterm 
d\ + d^. In fact, the loop contributions cannot be separated from the short-range pieces in a 
model-independent way. However, one does not expect cancellations between the non-analytic 
dependence in the pion mass and the analytic dependence of the short-range part. Taking the 
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long-range part in Eq. (182) at [i = mjy as an estimate for the neutron EDM, and using relation 



( [122D with 5m N ~ 2.3 MeV [60J, 

(- 4 ) I ™ I A M) 



The expected scaling of A^ R 4) in Eq. ([79]) with M QC d ~ 27rF T gives |d n | ~ 0.8 GeV e|£|/M : 



2 



which is about two thirds of the direct NDA estimate for d\ and d±, Eq. (164). We therefore 
expect the latter to provide a better estimate for the neutron EDM, and 

e\£\ 

K|~1.2GeV^. (189) 

Note that the neutron EDM was computed from general PJT four-quark operators in Ref. 
[38] . where it was found to be \d n \ ~ 0.15 GeV e\C$\ in terms of C4, the coupling constant of 
the four-quark operators. We compare to our result by writing C4 = (4ir) 2 £/Mp, and find that 
it is about an order of magnitude larger than our estimate. This discrepancy can be traced 
[63] to the use of a relativistic formalism in Ref. [3S]. In this reference, a large contribution 
to the nucleon EDM comes from pion loops in which the photon couples to the nucleon via 
the nucleon anomalous magnetic moment, while .PJT is provided by the coupling g\. While a 
chiral logarithm log m 2 / m 2 N appears with a suppression of m 2 /m 2 N , the contributions from loop 
momenta \k\ ~ give rise to a large EDM, which does not have such a suppression. In the 
heavy-baryon formalism, the nucleon magnetic moment does contribute to the nucleon EDM, 
but only at NNLO, and should reproduce the logm^ dependence of the relativistic calculation 
[63] . Unfortunately, xPT cannot be trusted for momenta where the nucleon is relativistic. Thus, 
we interpret the result of Ref. [48J as a model-dependent estimate of the size of the counterterms 
di, d\ and do, which, in this calculation, appears to be somewhat larger than NDA. 

10 Nuclear effects 

So far we have focused on interactions involving at most one nucleon. Power counting is more 
complicated for processes involving more than one nucleon: infrared enhancements due to purely 
nucleonic states require a resummation of leading interactions [23J , and this resummation in turn 
leads to the enhancement of certain CP-even operators that one would naively think are of high 



order [64J. Thus, the chiral index A in Eq. (62) is not particularly useful due to non-perturbative 
renormalization. 

On the other hand, ff operators should be treated perturbatively and are likely unaffected 
by this subtlety. It is convenient then to continue to organize them according to the index Ag 



in Eq. (73), where / can now take even values larger than 2. Since the index increases with 
increasing /, the operators most likely to be relevant have / = 4. The most important four- 
nucleon operators have been constructed in Ref. [35J . We give here some of the details for 
dimension-six sources. The consequences for the -PJT nuclear potential induced by the qCEDM, 
XISs and qEDM can be found in Ref. [35] . Ref. [35] did not discuss the f*p potential stemming 
from the FQLR operator. In this section, after constructing the short-range NN and NNw 
operators induced by each PJT source, we remedy this omission. 
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10.1 xISs 

While in the pion-nucleon sector one needs at least either two derivatives or one insertion of the 
quark mass to generate operators from the xISs, in the nucleon-nucleon sector one derivative is 
enough, and it need not include a pion field: 

C [ ~jl 4 = CiNNd^NS^N) + C 2 NtN ■ V^NS^tN), (190) 

with scalings 

c^o(^). (1.1) 

We see that for the xISs, the chiral index of the LO / = 4 LECs is the same as that of the LO 
/ = 2 LECs, in contrast with sources that are chiral breaking. Thus, whereas for 9, qCEDM, 
and FQLR the potential is dominated by pion exchange (and photon exchange for qEDM), the 
potential for the xISs gets contributions of the same order from pion exchange and short-range 
NN interactions |35j. 

The first terms linear in the pion field appear one order higher, 

TTi 

A =4 = • (GiNS' x tNV u (NS u N) + G 2 NS^NV v {NtS u N)) , (192) 

with scalings 

a " = °(^)- < 193 > 

When the pion is attached to another nucleon, a three-nucleon force results. 

The interactions above are not affected by tadpole removal because the operators are chiral 
invariant. At NNLO we find additional -fjT nuclear interactions, which we do not construct. 
At this order the first contributions to the isovector ffi nucleon-nucleon interactions appear as 
well, originating from an insertion of the quark mass difference. 

10.2 qCEDM 

Because it is chiral breaking, the qCEDM generates pion interactions at lower order than short- 
range, purely nucleonic interactions. At A6 = 0, 

4°/=4 = --^tt- (%NtNNN + ^ 2 NrS fl NNS fl N) 



7T3 



F V D 

where, after tadpole extermination, the LECs scale as 



j 3 NN NN + ^NS^N NS^N) , (194) 



1»=0 + 73,.=0^^j. (195) 

One order down, we find the first short-range interactions contributing directly to nucleon- 
nucleon scattering, 

4j=4 = (l " |J) [CiNN d^NS»N) + C 2 NtN • V^NS^tN)] 

+ (fa - [C 3 NnNdp (NS"N) + C^NNV^ (NnSW)] , (196) 
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with the scaling 

Here, as before, the LECs of the operators induced by the isoscalar qCEDM get a contribution 
from the isovector qCEDM as well, proportional to e. At this order, for all sources, there appear 
operators that start with one or more pions. We do not list them. 

A comparison between the pion-nucleon and nuclear sectors shows that for qCEDM the most 
important fp pion-nucleon interactions are larger by a factor Mq CD /Q 2 than the short-range 
nucleon-nucleon interactions, implying that for these sources the nuclear potential is 
dominated by pion exchange [35J. This observation justifies, a posteriori, the assumption often 
made in the literature, that nuclear observables can be calculated solely from fp pion 
exchange. However, as we have seen, for xISs this assumption is not valid. 



10.3 FQLR 

As we saw already in the pion and pion-nucleon sectors, the isovector qCEDM and FQLR 
generate very similar interactions, if one neglects the more complicated pion structure of the 
FQLR. 

Here again, 

4r 2 }=4 = 4/=4 + ' • ' • ( 198 ) 

where ". . ." represent interactions with more pion fields, which are of little importance. The 
LECs scale as 



after tadpole removal. Similarly, 



4rj=4 = 43=* + • • • ' ( 20 °) 



with 

As for qCEDM, one-pion exchange is more important than these short-range interactions [35], 
which we do not pursue further. 



10.4 qEDM 

The qEDM, while sharing with qCEDM and FQLR the property that interactions involving at 
least one pion dominate |35j, produces a richer isospin structure for the short-range nucleon- 
nucleon interactions, due to the integration of hard photon. 
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At Ag = 3 we find interactions that transform as the tensor product (— d^Vi + d 3 W 3 ) 
e(^/6 + r 3 V2), 



45=4 = --^—tv-^NtNNN + ^NtS^NNS^N) 
; ' (73 AW NN + jaNS^N NS^N) 



F n D 
~F V D 



5 3i + (vrgvri - n 2 5 3i ) [j 5 NnNNN + ^NnS^NNS^N] , (202) 



with 



2 

Q! em 771,,. 



(203) 



One order higher we find 

45=4 = (l - J^) [C 1 NNd ll (NS' i N) + C 2 iVYiV • V^NS^tN)] 

+ (s 3i - [C 3 Nr t Nd^ (NS^N) + C A NNV^ (Nt^N)] 

+ I (l - J^Q [CsN-nNdu {NS^N) + CeNNV^ (N n S»N)] 

+C 7 (fa - Nr k NV, (NS^N) } [s a + ^JL - 

+ . . . , (204) 
where we have ignored operators that start with one or more pions, and 



4vr F%Mj,Mqcd J ' ' ^ 4vr F%M 2 M QC D 



a ™ 2 



Since these operators are all suppressed by a em /Air, their phenomenological impact is minimal, 
and we do not construct operators with higher chiral index. 

10.5 The tff potential from the FQLR 

In Ref. [35J it was found that for the qCEDM the two-body Tf^ potential is dominated by 
one-pion exchange (OPE), where one of the pion couplings is a/^T coupling from Sec. [6j Short- 
range nucleon-nucleon interactions play a role at NNLO. On the other hand, for xISs short-range 
interactions appear at LO, while for the qEDM both pion and photon exchange contribute to 
the LO potential, which, however, does not usually play an important role in the calculation 
of EDMs. For qCEDM, qEDM and xISs, few-body forces are suppressed. For FQLR, we 
will see that the forces are somewhat different. We will show that, in LO, the two-nucleon 
potential is dominated by OPE and it is very similar to the potential induced by the qCEDM. 
However, the first loop correction appears already at NLO. The most striking feature of the 
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Figure 3: Contribution from the purely mesonic -fjT interaction to the three-nucleon potential. 
Notation as in Fig. [2} 



FQLR is, however, that the three-pion vertex A^ 4 ^ induces a three-nucleon force at LO, 
which impacts the calculation of the EDMs of 3 He and 3 H. 

We write the potential in terms of the spin, isospin, and incoming (outgoing) momentum of 
nucleon i, rW, and pi (p!) respectively. We denote the transferred momenta by q} = pi —p[. 
In the two-nucleon case, q\ = —q\, while for three nucleons, q\ = — {q\ + <fe). The coordinate- 
space version of the potential can be obtained straightforwardly following the procedure of, for 
example, Ref. [35] . 

At LO in xPT, the two-nucleon potential is given by OPE, with T violation provided by 



the pion-nucleon couplings of Sec. 6.3 The LO two-nucleon potential is thus identical to that 



generated by the qCEDM, and it is given by 

t/(-3) ) _ -9a qi r (i) (2 ) _ -.( 2 )\ 



(51 + 9l, 



206) 



Since the couplings go and g[ are not independent, this potential depends on two low-energy 



constants, which we can choose to be g\ and g[ using Eq. (123). The isoscalar piece of the 
potential is suppressed by the smallness of the ratio of the nucleon mass difference and the 
sigma term. Unless g\ is unnaturally tuned to cancel the contribution of g[, the two-nucleon 
potential is mainly isovector. 



In Sec. 7.4 we showed that the isovector PNFF F2 receives loop corrections, and in particular 
non-analytic momentum dependence, at NLO. These corrections are proportional to the three- 
pion vertex Al R 4 \ which is related to g[ by Eq. (122), and are enhanced by a factor of ir 



with respect to NDA. Their contribution to the isovector potential can thus be sizable, and, 
furthermore, it comes with a different momentum dependence with respect to the LO OPE. 

The deuteron EDM and MQM induced by the FQLR at LO can therefore be read off the 
calculations for the qCEDM of Refs. [361 GI], after one accounts for the suppression of go. At 
NLO, the deuteron EDM and MQM receive contributions from the three-pion vertex A^^\ 
and the NLO corrections could be important. While we postpone a detailed analysis of the 
A^ R 4 ^ correction to the deuteron EDM, we expect the main qualitative conclusion of Ref. [36] 
to continue to hold — that is, that the deuteron EDM induced by an isospin-breaking source 
is expected to be significantly larger than the isoscalar nucleon EDM. 

The three -pion vertex A^ R ^ has even more striking consequences in the three-nucleon system. 
This interaction induces a three-nucleon potential, depicted in Fig. [3j which in the power 
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counting of Ref. [36J enters at the same order as the LO two-nucleon potential. The three- 
nucleon potential is of the form 

Vtvlni^) = ^£(rf^).r(3) +T f T (D. T (3) +T (3) T (i). T (2) 



■ ^ ff(2) , ^ g(3) , ^ + ^ 

(<ff + m£)(gf + ml)[(qi + <f 2 ) 2 + m|] 



(207) 



The size of the contributions from Eq. (207) can be compared to those from the LO two-nucleon 



potential (206) by using the power-counting rules outlined in Ref. [36J. These rules indicate that 



the three-body contributions scale exactly as the LO two-body contributions. Equation (207) 
provides a correction to the helion and triton EDMs, which, by power counting, is of the same 
size as the two-body terms calculated in Ref. |36j . While we expect that the main qualitative 
conclusions of Ref. [36] — that for an isospin-breaking source the EDMs of helion and triton are 
significantly different from the EDMs of their constituents — to hold, we leave a quantitative 
analysis to future work. 

Notice that, as always, power counting provides only a guide for what should be included in 
a calculation. In particular, the power counting for few-nucleon forces has not been significantly 
probed. The one adopted in Ref. [35], which is in line with the current understanding of the 
two-body system, agrees with the one in Ref. [65]. If, instead, the power counting of Refs. 



[23], [66] is employed, then Eq. (207) becomes instead an NLO correction. 



The discussion of the two- and three- nucleon potentials shows that for observables in light 
nuclei the three-pion vertex A^ R 4 \ which survives tadpole extermination because of the tensor 
nature of the FQLR, plays an important role. Another handle to disentangle the qCEDM and the 
FQLR could be provided by the measurement of the T-odd correlation coefficient in f3 decay, D 



|22| . Indeed, the Hi operator in Eq. (24), which generates at tree level the FQLR operator, also 
contributes to an operator that couples right-handed quarks to left-handed electron and neutrino, 
with coefficient of O(Si). The qCEDM contributions to the same operator are suppressed by 
Gf due to the need for the exchange of a W boson and, because of an extra chirality flip, by 
an insertion of the light-quark mass. Therefore, the contribution of the qCEDM to the operator 
scales as dQ^Gpfhl^K ~ Gpffv^So^/M^,, which is suppressed by a factor m 2 /Myy- The size of 
the D coefficient with respect to nuclear EDMs is therefore very different for isovector qCEDM 
and FQLR operators. In case that the observation of nucleon and deuteron EDM point to an 
isospin-breaking source, the measurement of D could help to discriminate between these two 
operators. 



11 Discussion and conclusion 

As we have seen, the chiral structure of the various sources is influential in the form and 
expected magnitude of the low-energy ff interactions. Even without going into detailed results 
for _|°jr hadronic observables we can draw some qualitative conclusions by looking at the La- 
grangian that we constructed. In Ref. [32] it was found that for the Standard Model 9 term, all 
LECs are proportional to negative powers of the scale Mqcd- The reasons for this are twofold: 
the 9 term (i) can be seen as a complex quark mass term, which brings in at least one power 
of m 2 , and (ii) transforms, once vacuum stability is imposed, as the fourth component of an 
50(4) vector, which means that a pion tadpole can be eliminated. Just like isospin violation 
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|54j . time reversal is therefore an accidental symmetry in the SM, in the sense that it would be 
somewhat suppressed (by at least one power of m n / 'Mqcd) even if 6 were not small. 

Extending the analysis to the dimension-six operators, we see that positive powers of Mqcd 
do appear, but they are of course overcompensated by two negative powers of the much larger 
Mt- The main result of tadpole extermination is that the vacuum misalignment signaled by the 
pion tadpoles causes the isoscalar operators, like go or Ci,2, to receive additional contributions 
from isospin-breaking sources, like the isovector qCEDM and the FQLR operator, of the same 
importance as the contributions of isoscalar sources. These contributions can be schematically 
obtained by replacing 5q before tadpole removal with 6q + e5^ + e£,MQ CD /m% in the power 
counting estimates of isoscalar operators. Shifts in isovector quantities are inconsequential in 
the sense that they do not change the dependence of the LECs on £, <5o,3, and w. It is in the 
case of FQLR that tadpole rotation is most important because, first, FQLR does not have an 
isoscalar piece and, second, a relatively important three-pion vertex survives. 

We notice that the different chiral properties and field content of the 9 term, qCEDM, and 
FQLR on the one side, and qEDM and xISs on the other, imply very different relations be- 
tween long-distance and short-distance effects. The 9 term, qCEDM, and FQLR all violate 
chiral symmetry and thus generate pion-nucleon interactions in which the pion couples to 
the nucleon non-derivatively. As a consequence, the first f^f pion-nucleon couplings appear in 
the Lagrangian two orders before short-range contributions to the nucleon EDM. For the nu- 
cleon EDM and EDFF, this fact implies that even though pion-nucleon couplings can only 
contribute to the nucleon EDM via loops, which bring in a m^/(27T-F vr ) 2 suppression, they are 
still as important as the short-distance operators [29JH3]- For light nuclei, such as the deuteron 
\o'l\ 154"] and helion |67[ 136] , the most important contribution to f^P electromagnetic moments 
comes from the ffi OPE potential — which causes the nucleus wavefunction to mix with states 
of different parity — unless the admixed component has quantum numbers that cause the dipole 
matrix element to vanish. The application of chiral EFT to study effects of 9 term and qCEDM 
in systems with A > 2 nucleons is thus particularly promising, since observables are likely to 
depend on few LECs from the / = 2 Lagrangian. Once these constants are fixed in experi- 
ments, one is in the position to make testable, model-independent predictions. For the FQLR, 
there appear additional complications for EDMs of nuclei due to a three-pion interaction, 
which induces an NLO correction to the two-nucleon potential and, more importantly, a LO 
three-nucleon potential. The contributions of this three-nucleon force to light-nuclear EDMs 
have, so far, not been calculated. 

For the chiral-invariant sources, instead, the pion-nucleon coupling appear in the La- 
grangian at the same order as short-distance nucleon EDM operators. This happens because it is 
not possible to write a chiral-invariant pion-nucleon coupling with only one derivative. The 
first chiral-invariant pion-nucleon coupling must have two derivatives, while nonderivative 
couplings can be generated by considering the combined effects of chiral-invariant ff sources 
and the chiral-breaking quark mass; in any case, pion-nucleon couplings receive a further sup- 
pression of Q 2 /Mq CD . It should be noted that this difficulty does not affect the nucleon-photon 
and nucleon-nucleon sectors, where operators with a minimal number of derivatives can be 
constructed. The consequence for the nucleon EDM is that it is dominated in this case by short- 
distance contributions. For light nuclei, ff corrections to the wavefunction now are not only 
due to ff pion exchange, but also to short-distance nucleon-nucleon interactions. In general, 
the increased role of short-distance interactions in the case of xISs reduces the predictive power 
of our analysis, because of the appearance of more LECs. 
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Finally, pion physics is suppressed also in the case of -fjT from the qEDM. In this case 
the suppression comes from the need to integrate out the photon to produce purely hadronic 
operators, which leads to a factor of a e m/4vr. In this case, EDMs of light nuclei are dominated 
by the nucleon EDM. 

Traditionally, when discussing parity and time-reversal violation between nucleons and pions, 
three non-derivative pion-nucleon interactions are considered on the same footing: 

£?>tv = ~^-Nt ■ 77 N - ^-7T 3 NN - ^7T 3 Nt 3 N. (208) 
When one takes into account the chiral properties of the fundamental sources of fp, the 9 



term and the dimension-six sources in Eq. (54), this picture changes. First, these interactions 



are accompanied by others with further pion fields. Second, and more important, the three 
interactions are not of the same size. For example, in the case of the chiral-symmetry breaking 
but isospin-conserving 6 term, at LO only go appears [3]. The couplings g\ and g^ are respectively 
suppressed by two and three powers of m^/MQCD [22], although, once values from the connection 
to isospin violation in the quark masses [32] are taken into account, go is numerically smaller than 
expected [38J . If one wants to study observables sensitive to g\ one needs to take into account the 
full NNLO Lagrangian constructed above, which includes, apart from gi, also derivative pion- 
nucleon and multi-pion-nucleon interactions. An example of this is the fp two-nucleon potential 
which, at the order where g\ appears, has a rich and non-trivial momentum dependence in the 
isoscalar channel [35j . 

For qCEDM, arguably the most natural case is the one where the isoscalar and isovector 
components are of similar size, |<5o| — |^3|- In. this scenario the go and gi interactions appear 
at the same order. The third pion-nucleon coupling gi comes in two orders down in the chiral 
expansion. The case of a dominant isoscalar qCEDM, |<5o| 3> | <5^3 1 , generates an identical low- 
energy Lagrangian as the term, making it impossible to separate these two scenarios from 
low-energy fp observables alone. To this goal, more input from techniques like lattice QCD is 
needed. The appearance of a dominant isovector qCEDM, l^j 3> |<5q|, implies that go and g\ 
are approximately of the same order, although the former is expected to be somewhat smaller 
due to the extra factor e. 

The pattern of non-derivative pion-nucleon interactions in the case of the chiral- and isospin- 
symmetry-breaking four-quark operator FQLR is very similar to that of a dominant isovector 
qCEDM. The first difference is the appearance of interactions involving multiple pions, which are 
hard to isolate. A second difference is the appearance of fp mesonic operators which, through 
loops, add to the ffi PNFF momentum dependence, however at subleading order. The same 
mesonic operator induces a LO three-nucleon potential, which could be relevant for nuclei with 
A > 3. Regardless, separating the FQLR from an isovector qCEDM with hadronic observables 
requires very precise measurements and is, at this point, not very likely. The search for T 
violation in (5 decay could provide additional clues |22j . 

The xISs — gCEDM and two four-quark interactions — give rise to a similar hierarchy between 
the non-derivative pion-nucleon couplings as an isovector qCEDM, although in this case it is g\ 
that is expected to be smaller than go by a factor e. The LEC gi is suppressed by only one power 
of mjr/MQCD- Also, at the same order as go and g%, derivative pion-nucleon interactions appear. 
Within our framework a separation of the different chiral-invariant operators themselves is not 
possible. For that more advanced techniques than NDA are required to estimate the size of the 
LECs. 
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Table 1: The LO scaling of important LECs for the different sources. The g~i are the non- 
derivative 7riV couplings, do and d\ isoscalar and isovector short-range nucleon EDMs, cq and 
ci isoscalar and isovector magnetic "firN interactions, and C\ 2 and isoscalar and isovector 
AW interactions. (The + in some of the entries should not be taken literally but only as an 
indication that the LECs get contributions from two sources.) 



Source 


6 qCEDM FQLR qEDM xISs 


<?o 
9l/9o 
92/90 


77 m " /7- 7- \ 772^ Afnnn JVI^ nr , , , ^ TnP Adrtcrt ITi" Mnnn 
r "4 h 1 1 

A^qcd e 

rv 771™ Ao rv ^CICD 
«em c u o " U 7T f~ 7r u o "em ^ 


do 

d\/do 
co /do 


r/9 m * r(Sr 1 £ a l m - r c M QCD r m 2 r7« MqCD 
m qcd 1 ° ° 3) M^M QCU C ^ M 2 C " U M|M QCD CW M| 

1 1 1 | 1 

111 1 m * 

M 2 

JU QCD 

111 £3 "4 


Cl,2 
(73,4/(71,2 


"4 /X 1 ^-X \ m i £ M QCD fx 1 5 \ttcm "4 „„ M QCD 
^i^M^ (^0 + ^3) F 2 M , Mqcd tpT^ (So + S 3 ) 47r F 2 M | Mqcd ™^ 

r "4 <5 3 1 1 m 2 



Only for qEDM are the three non-derivative pion-nucleon couplings expected to be of the 
same order. But in this case all these couplings are pretty unimportant. The scalings of these 
couplings for all sources are summarized in the first three rows of Table [TJ 

The fact that different sources of f^p are responsible for different hierarchies between the 
non-derivative fp couplings has important implications for the -fjT electromagnetic moments of 
nuclei, of which the best example is the deuteron EDM. In calculations where the non-derivative 
pion-nucleon interactions are assumed to be equally sized, it is found that the deuteron EDM 
is dominated by g\ and is larger than the nucleon EDM [62J. The MQM gets contributions of 
similar size from go and g±. Similarly, it is found that for the helion EDM both go and gi are 
important [67]. Taking into account the chiral properties of the fundamental ff sources, these 
conclusions only hold for the isovector qCEDM and FQLR, since for these sources go and g\ 
interactions appear at leading order. More generally, we see no particular reason to insist on 
including 52 for any source, at least in calculations of light systems. 

Instead, what our analysis shows is that one should include other types of operators, 
such as nucleon-photon, pion-nucleon-photon, and nucleon-nucleon interactions. In Ref. [36], 
we argued that the EDMs of light nuclei depend on six different LECs. In addition to the two 
ff pion-nucleon interactions go and gi, they consist of the short-range isoscalar and isovector 
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nucleoli EDMs do and d\ , and the short-range isoscalar nucleoli- nucleon interactions Co and 
C\. In Ref. [36] the FQLR source was not considered. As we saw above, the surviving three-pion 
vertex introduces a dependence of observables on h^J^ 1 . However, A^ 4 ^ can be eliminated in 



favor of go and CP-even LECs via Eq. (122), so we can still consider LO nuclear EDMs to be 
expressed in terms of those six LECs. Higher electromagnetic moments, such as the deuteron 
MQM [l3l [37], can depend on additional interactions such as the isoscalar PJT pion-nucleon- 
photon vertex cq- We summarize the size of these important LECs (extended with some others) 
in Table [Q 

The Lagrangians derived in this work for each of the dimension-six PT sources serve as the 
basis for the calculations of hadronic and nuclear observables. In fact, the PT moments of the 
nucleon j29[ S3] [33] and of light nuclei [Ml ESI EH [38] have already been calculated with some of 
the interactions constructed here and, for the theta term, in Ref. |32|. Other observables such 



as the EDMs and Schiff moments of heavier nuclei can now be tackled with the same method. 
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